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FOREWORD 


This  section  presents  a detailed  account  of  the  project  activities  in  the  area  of 
hydraulic  noise.  Specific  test  procedures  are  recommended  for  determining  the  fluid- 
borne  noise  generation  potential  of  hydraulic  pumps  and  the  acoustical  performance 
characteristics  of  fluidborne  noise  attenuators.  Experimental  data  are  shown  to  indicate 
typical  results  that  are  obtained  with  the  proposed  test  codes. 


► 


TABLE  OF  CONTENTS 

Chapter 

Paqe 

1 

INTRODUCTION 

1-1 

1 1 

FLU  1 DBORNE  NOISE  ATTENUATOR  EVALUATION 

1-3 

T heory 

1-3 

Input  Impedance 

1-7 

Output  Impedance 

1-7 

Pressure  Drop 

1-8 

Frequency  Ranqe 

1-8 

T est  Procedure 

1-10 

T est  System 

1-12 

Test  Results 

1-1  3 

1 1 1 

PUMP  FLU  1 DBORNE  NOISE  GENERATION  POTENTIAL 

1-17 

Test  Procedure  Concept 

1-18 

T est  Systems 

1-26 

IV 

COMPUTER  PROGRAM 

1-37 

V 

CONCLUSIONS  AND  RECOMMENDATIONS 

1-47 

REFERENCES 

1-49 

Appendix 

A 

METHOD  FOR  EVALUATING  THE  PERFORMANCE  OF 
A FLUID 

1-51 

B 

PUMP  FLU  1 DBORNE  NC  1 SE  CHARACTER  1 ST  ICS 

1-57 

C 

PUMP  FLOW  RIPPLE  EVALUATION  TEST 

1-59 

D 

PUMP  IMPEDANCE  EVALUATION  TEST  PROCEDURES 

1-61 

E 

THE  OSU  SIX  TRANSDUCER  DATA  REDUCTION 
PROGRAM 

1-65 

LIST  OF  TABLES 


T abl  e 

Page 

2-1 

Table  of  Test  Frequencies  Using  Test  Plan  for  Obtaining 
15  Data  Points  for  Evaluating  Attenuator  Performance 

1-10 

3-1 

Slope  of  Flow  Ripple  Line  and  Statisfical  Information  for 
Two  Harmonics  and  Two  Pressures 

1-31 

3-2 

Reflection  Factor  and  Phase  Angle  for  the  Orifice 

1-31 

3-3 

Reflection  Factor  and  Impedance  of  Pump  OSU-NP-43 
for  Various  Frequencies 

1-33 

A-1 

Test  Conditions  Accuracy 

1-52 

A-2 

Data  Table  for  FBN  Attenuator  Test 

1-54 

D-l 

Data  Recording  Table  for  Pump  Impedance  Test 

1-63 

LIST  OF  FIGURES 

F i gure 

2-1 

Definition  of  Transmission  Loss  and  Relationship  Between 
Energy  Flows  to  and  from  a FBN  Attenuator 

1-6 

2-2 

Summary  of  Test  Procedure  Showing  Schematic  of  Measure- 
ment and  Loading  Circuit 

1-1  1 

2-3 

Transducer  Spacing  for  FBN  Attenuator  Performance 
Evaluation 

1-13 

2-4 

Transmission  Loss  Curve  for  FBN  Attenuator 

1-14 

2-5 

Reflection  Factor  Curve  for  a FBN  Attenuator 

1-14 

3-1 

Basic  Anechoic  Test  System  and  Equation  for  Evaluating 
the  Flow  Ripple  of  a Fluid  Power  Pump 

1-21 

3-2 

Basic  Test  System  Used  to  Evaluate  the  Source  Impedance 
of  the  Pump 

1-25 

3-3 

"Close-Coupling"  of  Pump,  Pressure  Transducer  and 
Orifice 

1-27 

3-4 

Reflection  Factor  for  Anechoic  Termination  Configuration 
Showing  the  Recommended  Maximum  Acceptable  Reflection 
Based  on  4dB  Standing  Wave 

1-28 

3-5 

Flow  Ripple  Versus  Pump  Speed  Showing  Experimental 
Values  and  Best  Fit  Line  Through  Origin  for  First  Harmonic 

1-30 

3-6 

Flow  Ripple  Versus  Pump  Speed  Showing  Experimental 
Values  and  Best  Fit  Through  Origin  for  the  Second  Harmonic 

1-30 

3-7 

Pump  Impedance  as  a Function  of  Frequency  (OSU-NP-43) 
Non-D  i mens  i ona  1 

1-35 

3-8 

Free  Volume  Flow  Ripple  as  a Function  of  Pump  Speed — 
First  Harmonic 

1-35 

3-9 

Free  Volume  Flow  Ripple  as  a Function  of  Pump  Speed — 
First  Harmonic 

1-36 

l-iii 


LIST  OF  FIGURES  (Continued) 


F i gure 


4-1 


4-2 


4-3 


A-1 

A-2 

A-3 


Schematic  of  Simple  Acoustic  System  Consisting  of  Source, 
Transmission  Line,  and  a Lead  Impedance,  Zj 

Six  Pressure  Transducers  Arbitrarily  Spaced  in  a FBN 
Measurement  System 

Simple  Flow  Diagram  of  Important  Computational  "Nested" 

Do-Loops  in  Program 

Pressure  Drop  Versus  Flow  Rate 

Transmission  Loss  Versus  Frequency 

Reflection  Factor  as  a Function  of  Frequency 


Page 


1-37 


1-41 

1-43 

1-55 

1-55 

1-55 


I — i v 


- FOREWORD  - 


This  report  was  prepared  by  the  staff  of  the  Fluid  Power  Research  Center  of  the 
Division  of  Engineering,  Technology  & Architecture,  Oklahoma  State  University  of 
Agriculture  and  Applied  Sciences.  The  study  was  initiated  by  the  Mobility  Equipment 
Research  and  Development  Command,  Fort  Belvoir,  Virginia.  Authorization  for  the 
study  reported  herein  was  granted  under  Contract  No.  DA AK02-75-C-01  37 . The  time 
period  covered  by  this  report  is  from  1 February  1976  to  31  January  1977. 

The  Contracting  Officer's  Representative  was  Mr.  Hansel  Y.  Smith,  and  Mr.  John  M. 
Karhnak  served  as  the  Contracting  Officer's  Technical  Representative.  In  addition, 

. Paul  Hopler  has  effectively  represented  the  Contracting  Officer  both  technically 
administratively  through  various  phases  of  this  contract.  The  active  participation 
Vssrs.  Smith,  Karhnak,  and  Hopler  during  critical  phases  of  work  contributed 
cantlv  to  the  overall  success  of  the  program. 

The  studies  represented  by  this  report  were  conducted  under  the  general  guidance 
of  Dr.  E.  C,  Fitch,  Program  Director.  The  details  of  each  study  are  presented  in  a 
self-contained  section  of  this  report.  The  titles  of  the  various  sections  together  with 
their  respective  Project  Managers  are  listed  below: 

SECTION  I . 

SECTION  I I . 

SECTION  III. 

SECTION  IV. 

SECTION  V. 

SECTION  VI  . 


HYDRAULIC  NOISE  ATTENUATION  — G.  E.  Maroney 
CYLINDER  STRUCTURAL  INTEGRITY  ASSESSMENT  — 
S.  K.  R.  Iyengar 

HYDRAULIC  SYSTEM  OPERATIONAL  SEVERITY 
ASSESSMENT  — R.  L.  Decker/S.  K.  R.  Iyengar 
ON-BOARD  MONITOR  STUDY  — R.  L.  Decker 
HYDRAULIC  SYSTEM  DIAGNOSTICS  — R.  K.  Tessman 
PUMP  CONTAMINANT  TOLERANCE  VERIFICATION  — 
L.  E.  Bensch 


CHAPTER  I 


I NTRODUC 1 ION 


The  need  to  control  the  noise  emitted  from  fluid  power  equipment  has  stimulated 
many  manufacturers  and  users  in  the  fluid  power  industry  to  take  a "fresh  look"  at 
noise.  There  has  been,  in  recent  years,  a more  intensive  look  at  the  sources  of  noise 
in  order  to  understand  and  model  these  sources.  Instead  of  concentrating  only  on  air- 
borne noise,  research  has  been  extended  into  the  areas  of  fluidborne  and  structure- 
borne  noise  since  the  control  of  these  two  forms  of  noise  can  have  a significant  impact 
on  the  overall  airborne  noise.  This  report  is  limited  to  the  area  of  fluidborne  noise 
(FBN). 

Fluidborne  noise  can  be  controlled  by  either  changing  the  source  of  the  noise 
(usually  a pump,  motor,  or  valve),  installing  a FBN  attenuator,  or  otherwise  modify- 
ing the  system.  However,  to  change  the  source,  it  is  necessary  to  understand  the 
important  parameters  that  constitute  the  fluidborne  noise  generation  potential  of  the 
source.  In  this  way,  say,  two  pumps  can  be  compared  to  see  which  one  has  the  largest 
(or  smallest)  fluidborne  noise  potential.  Also,  if  the  important  parameters  of  a noise 
source  are  determined  (either  analytically  or  experimentally),  then  it  is  possible  to 
predict  how  the  pump  will  affect  a given  system. 
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A similiar  procedure  can  be  used  for  FBN  attenuators.  The  use  of  attenuators 


in  fluid  power  systems  has  often  been  with  a "fly  by  the  seat-of-your-pants"  philosophy. 
The  person  installing  in  attenuator  in  a system  hopes  that  a given  attenuator  will  be 
the  "one"  that  will  work  for  his  system.  With  an  adequate  understanding  of  the  per- 
formance of  attenuators,  it  would  be  possible  for  manufacturers  and  users  to  compare 
attenuator  A with  attenuator  B and  predict  what  effect  they  would  have  on  the  system. 

Chapter  I I of  this  report  discusses  a fluidborne  noise  attenuator  evaluation  pro- 
cedure. The  basic  theory  and  important  parameters  are  outlined,  as  well  as  an  actual 
test  procedure.  This  procedure  is  applied  to  an  actual  attenuator  to  evaluate  its  per- 
formance characteristics  as  compared  to  the  theoretical  model. 

Chapter  I I I covers  the  area  of  FBN  sources,  in  particular  pumps.  A basic  test 
procedure  is  outlined  for  the  evaluation  of  the  FBN  potential  of  fluid  power  pumps  to 
include  source  impedance  and  flow  ripple.  These  are  compared  to  the  theoretical 
model . 

Chapter  IV  is  a brief  discussion  of  the  OSU  Six  Transducer  Data  Reduction  Com- 
puter Program.  This  program  has  added  tremendously  to  the  OSU-FBN  program  by 
providing  a more  accurate  estimation  of  standing  wave  parameters  such  as  incident 
pressure,  reflection  factor,  and  the  position  of  maximum  pressure. 
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CHAPTER  I I 


FLU  I DBORNE  NOISE  ATTENUATOR  EVALUATION 


Accompanying  the  increased  availability  and  use  of  fluidborne  noise  attenuators, 
there  is  an  even  greater  need  for  an  objective,  reliable  and  reputable  method  of  eval- 
uating the  effectiveness  of  these  attenuators.  Since  attenuators  are  designed  to  reduce 
the  fluidborne  noise  (FBN)  or  pressure  ripple  in  a system,  then,  to  critically  evaluate 
an  attenuator,  a fundamental  understanding  of  FBN  is  essential.  This  section  of  the 
report  will  outline  some  of  the  basic  equations  and  assumptions  needed  for  evaluating 
attenuators  outline  a test  procedure  for  FBN  attenuator  evaluation  and  present  test 
results  from  a sample  attenuator. 

THEORY 

| 

I 

, To  model  fluidborne  noise  or  wave  propagation  in  fluid  power  systems,  plane 

wave  theory  can  be  utilized  by  making  the  following  assumptions: 

1 . Wave  propagation  is  one-dimensional . This  condition  is  satisfied  as 
long  as  the  diameters  of  pipes,  conduits,  etc.  , are  less  than  one- 

, fourth  the  wave  length  of  the  highest  frequency  of  interest.  For 

instance,  in  a 2.54  cm  (1  inch)  i.  d.  pipe  containing  a fluid  having 
a sonic  velocity  of  1220  m/s  (4000  ft/sec),  plane  wave  theory  can  be 
used  up  to  frequencies  of  12,000  Hz. 
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2.  Viscous  Dissipation  is  negligible.  For  very  high  flow  rates  in  small 
diameter  tubes,  plane  wave  theory  has  to  be  modified  to  include  a 
viscous  dissipation  term. 

3 . The  amplitude  of  fluctuation  in  flow  and  pressure  is  small  in  com- 
parison with  mean  flow  and  pressure.  This  assumption  linearizes 
the  differential  equations  used  to  derive  plane  wave  theory.  If  this 
condition  is  not  satisfied,  as  it  may  be  in  some  fluid  power  systems, 
then  a nonlinear  model  must  be  used. 

Given  these  assumptions,  the  acoustic  power  flow  per  unit  of  area  in  an  infinitely 
long  pipe  is  the  rms  pressure,  p,  times  the  rms  particle  velocity.  The  rms  particle 
velocity  is  p/Z,  where  Z is  the  characteristic  impedance  of  the  pipe  which  is  defined 
as  pC/S,  with  S being  the  cross-sectional  area  of  the  pipe.  The  power  in  the  pipe 
is  [7]: 

w = S P2/ PC  (2-1) 

Infinitely  long  lines  seldom  exist  in  the  field  hydraulic  systems,  so  there  is  usually 

2 2 

a right  traveling  wave,  (S  p.  / pC),  and  a left  traveling  wave,  (S  p^  / pC).  The  net 
power  flow  is: 

W , - (S/PC)  (p.2  - P 2)  (2-2) 

net  i r 

where  p.  and  p^ , respectively,  are  associated  with  the  incident  and  reflected  pressure 
waves.  Eq.  (2-2)  can  be  shown  to  be  equivalent  to: 
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(2-5) 


net 


(S/  PC)  (p  ) (p  . ) 
max  min 


where  p is  the  maximum  pressure  (rms)  in  the  line  and  p is  the  minimum  pressure 
max  mm 

in  the  line.  The  ratio  p /p  is  defined  as  the  standinn  wave  ratio  (SWR).  The 

max  min 

followinq  equations  show  relationships  between  SWR , p.  , p , p , and  p . : 

J i r max  nun 


SWR 


P /p  ■ 
max  min 


(p.  + Pr)/(Pj  - Pr) 


(2-4) 


or  rearranging: 


pr/pi 


(SWR  - 1 )/(SWR  + 1) 


(p  - P - )/(p  + p . ) 

max  min  max  min 


(2-5) 


IMPORTANT  TEST  PARAMETERS 

A test  code  for  the  evaluation  of  fluidborne  noise  attenuators  should  consider  the 
following  characteristics:  transmission  loss,  input  impedance,  output  impedance,  and 
flow  resistance.  The  transmission  loss  indicates  how  much  of  the  incident  pressure 
ripple  is  emitted  at  the  outlet  of  the  attenuator.  The  input  impedance  indicates  how  much 
of  the  incident  pressure  is  reflected  by  the  attenuator  toward  the  source.  The  output 
impedance  indicates  how  much  of  the  incident  pressure  is  reflected  by  the  attenuator 
toward  the  source.  The  output  impedance  indicates  what  percentage  of  the  pressure 
incident  on  the  attenuator  outlet  is  reflected  back  to  the  downstream  portion  of  the 
system.  The  flow  resistance  of  the  attenuator  indicates  the  power  dissipated  by  a given 
value  of  mean  flow  through  the  unit.  These  four  characteristics  adequately  define  the 
performance  characteristics  needed  to  acoustically  rate  the  attenuator  and  to  compare 
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tiie  power  dissipated  by  the  attenuator  relative  to  the  acoustical  effectiveness  of  the 
unit . 


TRANSMISSION  LOSS 


Wt-^(P2-Pf) 
P,  


J U 

fluidborne 

NOISE 

ATTENUATOR 


w2*4?(p?- 

p„ 


p,20 


TL  = 10  log,0$jMd8) 

W,=-^-(pf) 

W,  =^-(pf) 


Fiq.  2-1  . Definition  of  T ransmission 
Loss  and  Relationship  Between  Energy 
Flows  to  and  from  a FBN  Attenuator. 


Figure  2.1  shows  the  relationship 

between  p.  and  p in  a section  of  conduit 
i r 

connected  to  a fluidborne  noise  attenu- 
ator. Extending  the  same  concept  to 
the  pipe  downstream  of  the  attenuator 
gives  a relationship  for  the  transmitted 
pressure,  p{ , and  the  reflected  pressure, 
p(r,  in  the  downstream  pipe.  One  per- 
formance parameter  for  the  attenuator 
is  transmission  loss: 


TL  = 10  log1()  (W./Wt)  (2-6) 

If  the  areas  of  the  upstream  and  downstream  conduits  are  equal  and  the  system  has 
equal  densities  and  sonic  velocities  in  both  sections,  then  it  follows  that: 


TL  = 10  log10  (S1p.2/p1C1)(P2C2/S2Pt2) 


TL  = 10  log1Q  (p^/p^) 


(2-7) 
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TL  = 20  log,  (p./p  ) 

J1o  i t 

Thus,  the  transmission  loss  can  be  determined  by  evaluation  of  experimental  data  using 
Eqs.  (2-6)  and  (2-7). 


INPUT  IMPEDANCE 


Ignoring  phase  relationships  for  this  report,  we  can  consider  only  the  real  part 
of  the  input  impedance  of  the  attenuator  that  is  the  absolute  value  reflection  coefficient, 
the  reflection  factor,  R.  The  reflection  factor  is  p^/p..  R can  be  evaluated  using 
Eq.(2-5).  The  data  for  evaluating  the  input  reflection  factor  is  the  same  data  used 
to  evaluate  the  transmission  loss. 


OUTPUT  IMPEDANCE 

If  the  attenuator  will  accept  flow  through  the  unit  from  either  direction,  it  is 
possible  to  evaluate  the  output  impedance  using  the  output  reflection  factor,  R . The 
data  for  calculating  Rq  is  obtained  by  inverting  the  attenuator  in  the  test  circuit  and 
obtaining  the  necessary  measurements  by  repeating  the  test  procedure  used  for  obtain- 
ing the  transmission  loss  and  the  imput  impedance. 

If  the  pressure  drop  through  the  attenuator  is  less  than  or  equal  to  an  equivalent 
length  of  tubing,  then  the  attenuator  is  essentially  non-dissipative.  For  non-dissipative 
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units , 


2 

p^"-,  and  Eq.  (1-7)  can  be  rewritten: 


TL  = 10  log.n  (p.2/(p.2  - p 2) 

I U i i r 

TL  = 10  log, Q (1/(1  - R2) 


(2-0) 


Thus,  the  reflection  factor  for  non-dissipative  elements  can  be  written  in  terms  of  the  TL 


R2  = 1 - 10-(TL/10) 


(2-9) 


PRESSURE  DROP 

The  efficiency  of  the  attenuator  can  be  evaluated  by  ratioing  the  attenuation  of  the 
unit  to  the  pressure  drop  of  the  unit  at  the  desired  flow.  To  obtain  an  estimate  of  the 
efficiency  of  the  attenuator,  data  are  needed  which  show  the  pressure  drop  versus  flow 
of  the  unit.  These  data  can  be  obtained  using  an  acceptable  test  procedure. 


FREQUENCY  RANGE 

One  reason  for  uainq  a fluidborne  noise  attenuator  is  to  reduce  the  amplitude  of 
system  pressure  ripple  at  the  pumping  frequency.  Since  most  conventional  pumps  have 
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between  9 and  15  pumping  elements  (n)  and  are  usually  operated  above  600  revolutions 
per  minute  (rpm)  (N)  but  at  or  below  2500  rpm,  a fundamental  pumping  frequency  ranqe 


f,  = Nn/60 

flow  = (600)  9/60  = 90  Hz  (2-9) 

fhigh  = (2500)15/60  = 625  Hz 

If  the  proposed  procedure  requires  that  the  first  three  pumping  harmonics  are  to  be 
measured,  this  will  accommodate  a frequency  range  of  interest  spanning  100  Hz  and 
2000  Hz.  Constraining  the  frequency  range  to  100  Hz  - 1000  Hz  allows  testing  to  be 
conducted  in  the  average  fluid  power  laboratory,  using  maximum  pump  speeds  in  the 
vicinity  of  2300  rpm  for  a pump  with  nine  pumping  elements,  thus  providing  information 
about  the  attenuator  over  a frequency  range  which  is  important  to  the  largest  number 
of  users . 

Because  the  pressure  level  at  the  fundamental  pumping  frequency  is  usually  more 
repeatable  than  the  levels  at  higher  harmonies  and  because  most  attenuators  will  prob- 
ably respond  best  to  the  highest  pressure  level  in  the  system,  it  is  recommended  that 
data  for  attenuator  evaluation  be  taken  at  the  first  three  pumping  frequencies.  Data 
at  individual  frequencies  are  obtained  by  varying  the  pump  speed.  At  a minimum,  the 
procedure  requires  measuring  attenuator  performance  at  the  eleven  third-octave  fre- 
quencies between  100  and  1000  Hz. 
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TEST  PROCEDURE 


F iqure  2-2  outlines  the  test  procedure.  (Appendix  A contains  a more  detailed 
and  specific  procedure.)  Table  2-1  gives  the  test  frequencies  that  should  be  used  to 
evaluate  the  FBN  attenuator  if  the  first  three  harmonics  of  a pump  are  used.  This 
yields  data  points  at  the  eleven  third-octave  frequencies  between  100  and  1000  Hz 
plus  four  overlapping  data  points. 


TABLE  2-1 

Table  of  Test  Frequencies  Using  Test  Plan  for  Obtaining  15  Data  Points  for  Evalu- 
ating Attenuator  Performance 


Desired 
1/3  Octave 
Center 

TEST  NUMBER 

F requency 
(Hz) 

1 

2 

3 

4 

5 

100 

100* 

125 

125* 

160 

160* 

200 

200 

250 

250 

(267)* 

315 

(300) 

(320) 

315* 

400 

(375) 

500 

(480) 

(533) 

630 

630 

800 

800 

1000 

(945) 

! 

• * Fundamental  pumping  frequencies 

NOTE:  Frequencies  in  parentheses  deviate  from  desired  frequencies. 
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Step  1 
Step  2 
Step  3 
Step  4 
Step  5 
Step  6 

Step  7 


Verify  anechoic  termination  and  measure  "tare1’  P. 
Measure  pressure  drop  across  attenuator. 

Determine  energy  incident  to  attenuator  (p.)  . 2 

Determine  energy  transmitted  through  attenuator  (p.)  . 
Measure  sonic  velocity. 

T.  L.  = 20  log1Q  P/Pt- 

Calculate  pressure  drop. 


Fig.  2-2.  Summary  of  Test  Procedure  Showing  Schematic  of  Measurement  and 
Loading  Circuit 
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TEST  SYSTEM 


The  critical  components  of  the  test  circuit  are  illustrated  in  Fig.  2-2.  The  max- 
imum flow  of  the  pump  at  a speed  to  provide  the  1000  Hz  test  pressure  ripple  should  be  - 
less  than  the  maximum  allowable  flow  of  the  attenuator.  The  hydraulic  lines  from  the 
pump  to  the  load  valve  should  be  steel  tubing  of  the  source  diameter,  connected  where 
necessary  with  unions  that  minimize  discontinuities  in  the  line.  The  transducer  sensing 
surfaces  should  be  tangent  to  the  inside  diameter  of  the  steel  conduit. 

The  range  of  attenuator  sizes  that  can  be  accommodated  by  a given  system  can 
extend  above  and  below  the  system  conduit  size  if  controls  are  placed  on  the  transi- 
tion lengths  of  porting  adapters.  Expansion  chamber  design  curves  [7]  indicate  that 
porting  adaptor  transitions  with  lengths  less  than  0.25  X should  not  affect  the  test  re- 
sults. A 

\\ 

A test  system,  such  as  shown  in  Fig.  2-2,  which  has  several  meters  of  conduit 
between  the  measurement  region  and  the  load  valve  will  provide  a quasi-anechoic  termin- 
ation. This  means  that  the  variations  of  measurements  as  a function  of  x will  be  greatly 
reduced.  I f enough  conduit  is  used  so  that  R "0,  then,  from  Eq.  (2-1),  p(x)  = p. . 

If  the  measurement  variations  are  less  than  2 dB  at  any  frequency,  the  reading  can  be 
averaged  on  a dB  basis  to  obtain  a value  representative  of  p..  Since  large  variations 
of  the  pressures  will  introduce  the  possibility  of  significant  errors  in  obtaining  p. , it 
is  reasonable  to  set  some  upper  limit  on  the  range  of  dB  readings  at  a given  operating 
condition.  The  test  system  should  meet  the  qualification  that,  at  any  test  condition, 
the  range  of  dB  readings  shall  not  exceed  4 dB.  This  limit  still  allows  an  appreciable 

standing  wave;  the  standing  wave  ratio  is  1 .6,  or,  p is  equal  to  1 .6  p . . 

J ^ max  min 


mar: 


The  transducer  configuration  should  include  at  least  three  pressure  transducers. 
It  is  recommended  that  six  transducers  be  used  with  the  spacing  shown  in  Fig.  3-3. 
The  configuration  shown  in  Fig.  2-3  provides  better  accuracy  than  three  transducers. 
The  spacing  is  based  on  the  Renard  Series. 


5 m-- 


0 25m-t~$-'V075m-h 


1 00m*t* 

— 5- 


2 50m  — 


“ft 


Fig.  2-3.  Transducer  Spacing  for  FBN  Attenuator 
Performance  Evaluation 


TEST  RESULTS 

The  attenuator  selected  for  evaluation  was  an  expansion  chamber,  which  is  a 
device  with  an  increase  in  diameter  at  the  inlet  and  a corresponding  reduction  in  dia- 
meter at  the  outlet.  The  unit  chosen  for  the  test  was  symmetrical  as  described  in 
Appendix  E.  The  unit  was  tested  for  attenuation  characteristics,  input  impedance 
and  flow  resistance.  The  results  of  the  tests  are  described  in  the  following  paragraphs. 
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REFLECTION  FACTOR 


I ransmission  Loss 


'RAMSMSSK*  LOSS  Of 
ATTf  NUATOft 


Fiq.  2-4.  Transmission  Loss  Curve 
for  F BN  Attenuator 


REFLECTION  FACTOR  FOR 
F8N  ATTENUATOR 


FREQUENCY 

Fig.  2-5.  Reflection  F actor  Curve  for 
a FBN  Attenuator 


Fig.  2-4  compares  the  results  of 
the  transmission  loss  evaluation  with 
the  theoretical  transmission  loss  for 
the  expansion  chamber.  The  experi- 
mental transmission  loss  was  evaluated 
by  using  the  OSU  Six  Transducer 
Data  Reduction  Program.  The  trans- 
ducer configuration  was  the  same  as 
suggested  in  Fig.  2-3. 

The  scatter  using  the  six  trans- 
ducers is  significantly  less  when  com- 
pared to  the  data  scatter  obtained  from 
three  transducers.  The  six  transducer 
technique  appears  to  give  more  accurate 
values  of  the  transmission  loss  for  the 
frequencies  measured.  The  procedure 
for  determining  the  transmission  loss  of 
a FBN  attenuator  is  given  in  Appendix  A. 

Fig.  2-5  shows  the  reflection  factor 
for  the  input  to  the  expansion  chamber. 
These  values  were  determined  using  the 
OSU  Six  T ransducer  Data  Reduction 
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Program.  For  a non-dissipative  F BN  attenuator,  the  reflection  factor  of  the  attenuator 
can  be  deduced  from  the  transmission  loss  with  the  use  of  Eq.  2-9.  These  values  are 
also  shown  in  Fig.  2-5. 

As  is  shown  in  Fig.  2-5,  it  appears  that  the  OSU  Six  Transducer  Data  Reduction 
Program  can  directly  provide  accurate  information  about  the  reflection  factor  up  to 
approximately  325  Flz.  Above  325  Flz  the  program  (and  data)  seem  to  give  values  of  R 
that  are  smaller  than  expected.  This  might  be  more  of  a limitation  on  the  measurement 
technique  itself  than  the  computer  program.  The  data  was  not  taken  in  "real  time"  so 
that  any  sonic  velocity  changes  during  the  tests  would  have  affected  the  measured  pressure 
values.  This  would  affect  the  fit  that  the  computer  program  would  have  for  the  data. 

Another  possible  source  of  error  is  the  decrease  in  wavelength  with  a correspond- 
ing increase  in  frequency.  For  instance,  at  100Hz  there  are  six  transducers  covering 
slightly  less  than  half  a wavelength.  At  1000Hz  there  are  five  wavelengths  between  the 
first  and  last  transducer.  With  5 "peaks"  and  "valleys"  it  is  harder  to  locate  the  exact 
location  of  maximas  and  minimas  since  they  are  so  much  closer  together.  The  ideal, 
of  course,  would  be  to  have  infinite  number  of  transducers  in  order  to  directly  measure 
the  maximas  and  minimas.  However,  this  is  not  realistic.  More  study  should  be  done 
to  determine  the  maximum  and  average  error  an  investigator  might  expect  from  using 
the  six  stationary  transducer  technique. 


1-15 


CHAPTER  I I I 


PUMP  FLU I DBORNE  NOISE  GENERATION  POTENTIAL 


With  the  growing  concern  for  the  reduction  of  the  overall  airborne  noise  level  in 
hydraulic  systems,  attention  has  begun  to  be  focussed  on  understanding  the  interaction 
that  contribute  to  the  overall  airborne  noise  level.  In  particular,  this  has  meant  a 
more  "in-depth"  look  at  fluidborne  noise  and  structureborne  noise  [1  , 3,  4,  5,  7,  9,  10]. 

A major  source  of  fluidborne  noise  in  hydraulic  systems  is  the  positive  displacement 
pump.  Since  the  fluidborne  noise  level  can  contribute  significantly  to  the  overall 
airborne  noise  level  in  a system,  the  need  has  developed  to  adequately  model  and  rate 
components  that  generate  fluidborne  noise;  namely,  the  positive  displacement  pump. 

The  fluidborne  noise  (FBN)  generated  by  fluid  power  pumps  has  been  the  subject 
of  several  papers  and  reports  in  recent  years  [2,  3,  4,  5,  9,  10J.  Ichikawa  and 
Yamaguchi  [2]  developed  a model  for  the  variation  in  flow  and  pressure  (FBN)  generated 
by  a gear  pump.  From  the  data  given,  the  model  proved  satisfactory  for  predicting  the 
flow  ripple  of  the  gear  pump.  They  consider  the  leakage  (or  source)  impedance  in 
calculations  to  determine  the  pressure  ripple  delivered  to  the  line. 

Wi  I likens  [3]  presents  a model  for  a gear  pump  and  presents  data  to  support  his 
model.  He,  also,  shows  that  the  pressure  ripple  measured  is  dependent  on  both  the 
flow  ripple  generated  by  the  pump  and  the  leakage  impedance  of  the  pump. 

' ' ' "* " 
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Unruh  [5,  p.2]  recognizes  the  need  to  ultimately  relate  the  measured  pressure 
ripple  back  to  the  flow  ripple  of  t fie  pump.  He  states:  "...  the  basic  phenomenon 
which  actually  needs  to  be  established  for  determining  the  fluidborne  noise  character- 
istics of  a pump  is  the  amount  of  flow  ripple  that  a pump  generates  due  to  its  inherent 
construction  characteristics  . . . ." 

If  a pump  can  be  considered  as  being  analogous  to  a constant  current  source  in 
an  electrical  circuit,  then  the  pump  will  deliver  a constant  flow  ripple  to  the  system  at 
a given  load  condition.  Just  as  the  load  on  a constant  current  source  determines  the 
voltage  in  the  circuit,  so  the  load  conditions  on  a fluid  power  pump  will  determine  the 
magnitude  of  the  pressure  ripple  in  the  system.  The  pump  also  has  a source  impedance 
analogous  to  that  of  a constant  current  source.  For  a pump,  the  source  impedance  will 
depend  on  the  geometry  (or  volume)  of  the  pump  casing  and  the  amount  of  leakage  flow  [2,3]. 

The  two  important  parameters  in  defining  an  electrical  current  source  are  the 
strength  of  current  and  the  source  internal  impedance.  Similarly,  for  a hydraulic 
pump,  the  inherent  flow  ripple  of  the  pump  and  the  source  impedance  are  needed  to  ade- 
quately describe  the  pump  in  order  to  model  and  predict  the  fluidborne  noise  level  in  a 
hydraulic  system.  A general  approach  for  evaluating  these  parameters  of  fluid  power 
pumps  is  discussed  in  the  next  section. 

TEST  PROCEDURE  CONCEPT 

To  determine  the  flow  ripple  and  impedance  of  a fluid  power  pump,  several  different 
measurements  must  be  taken.  The  flow  ripple  can  be  determined  by  measuring  the  pres- 
sure ripple. 
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The  basic  concept  behind  tine  proposed  test  procedure  for  determining  the  fluidborne 
noise  potential  of  fluid  power  pumps  involves  two  phases: 

1 . Determining  the  Flow  Ripple  - With  today's  technology,  it  is  not  industrially 
practical  to  directly  measure  the  flow  ripple  of  fluid  power  pumps.  Pressure 
ripple  measurements  must  be  taken  near  the  pump  with  an  anechoic  termination, 
and  then  converted  to  flow.  This  is  discussed  in  the  section,  F low  Ripple 
Determination . 

2.  Determining  the  Pump  Impedance  - This  phase  involves  coupling  a known  load 
impedance  within  1/20  of  a wavelength  of  the  pump  and  measuring  the  pressure 
ripple  between  the  pump  and  known  load  impedance.  From  this  information,  the 
source  impedance  can  be  calculated  as  is  discussed  in  the  section.  Pump 
Impedance  Determination.  There  are  two  distinct  steps  in  the  phase.  The 
first  is  the  evaluation  of  the  load  impedance  and  second,  the  evaluation  of  the 
source  impedance. 

Flow  Ripple  Determination 

As  stated  in  the  previous  section,  determination  of  the  flow  ripple  of  the  pump 
involves  taking  pressure  ripple  measurements  near  the  pump,  which  is  connected  to 
anechoic  termination.  An  anechoic  termination  would  be  a termination  which,  ideally, 
would  produce  no  reflections.  The  noise  in  the  fluid  would  propagate  down  the  tube  and 
not  be  reflected  back  toward  the  source.  Since  fluid  power  pumps  have  a wide  range  of 
pumping  freguencies,  an  anechoic  termination  needs  to  be  "non-reflective"  over  a wide 
range  of  freguencies. 


1-19 


The  difficulty  in  obtaining  a true,  non-reflecti ve  termination  has  been  discussed 
by  Heymann  [b].  He  states:  "It  is  . . . difficult  to  obtain  a practical  reflection-free 
termination  in  liquid  systems,  most  approaches  realize  at  best  a termination  in  which 
reflections  . . . are  weak  . . . . " 


One  approach  to  constructing  an  anechoic  termination  is  to  use  very  long  lengths 
of  line  so  that  the  pressure  ripple  is  eventually  dissipated  as  it  travels  down  the  line. 
For  pipe  diameters  of  one  inch  and  flows  above  3.5  gpm,  a hydraulic  circuit  line  length 
of  between  35  - 70  meters  is  required.  For  larger  diameters,  the  length  required 
will  be  longer,  since  the  frictional  dissipation  will  decrease.  Also,  for  larger  flows, 
the  length  of  line  required  will  decrease  as  the  flow  increases.  Generally,  as  the 
frequency  increases,  the  length  required  should  decrease.  The  dissipation  term  is 
thus: 


D = F(d,  Q,  f) 


(5-1) 


where:  D 

d 
Q 
f 


dissipation  term 
pipe  diameter 
flow 

frequency 


The  basic  concept  of  measuring  pump  pressure  ripple  in  a system  with  an  anechoic 
termination  is  shown  in  Fig.  3—1  . In  Fig.  3—1  , a pump  is  connected  to  an  anechoic 
termination  ahead  of  the  system  load  valve  and  filter.  If  there  are  no  reflections  at  the 
termination,  then  the  pressure  at  any  point  in  the  hydraulic  line  will  be  the  same,  since 
there  will  be  no  standing  waves.  In  such  a case,  the  impedance  at  any  point  in  the  line 
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P(x)  = P, 


Qo=  P/Z0 


Fig.  3-1.  Basic  Anechoic  Test 
System  &■  Equation  for  t valuatinq 
the  Flow  Ripple  of  a Fluid  Power 
Pump 


is  just  the  characteristic  impedance,  Z . Thus,  the  flow  ripple  in  an  anechoic  termina- 
tion is  related  to  the  pressure  as: 


where: 


°o  = V Zo  ($-2) 

P = measured  pressure  ripple 

o 

= characteristic  impedance  of  the  line  - pC/S 
p = density  of  the  fluid 

C = sonic  velocity 

S = cross-sectional  flow  area 


The  general  relationship  between  the  free  volume  flow  ripple 
pump,  Qj,  and  fl^w  ripple  measured  in  the  line,  0^  , is  given  by 


generated  by  the 
[18]: 


O. 


Z, 


Zs  + ZL 


(3-3) 
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where  Z and  Z^  are  the  impedances  of  the  pump  and  load,  respectively.  If  an  ane- 

choic  termination  is  used  as  the  load,  then  Z,  Z and  Q.  O , in  which  case: 

L o L o 


Q 

o 


"S 


+ Z 

o 


(3-4) 


There  is  a similar  relationship  between  the  pressure  ripple  measured  in  a line,  , 
and  the  pressure  ripple  generated,  P^,  by  the  pump.  This  is  [18]: 


P = P /L 
L B L 


(3-5) 


Z 


S + 


For  an  anechoic  termination,  this  becomes: 


O 


B 


(3-6) 


S + 


A look  at  Eq.  (3-4)  shows,  that  if  Z^>>  ZQ*  then  the  flow  ripple  measured  in  an 
anechoic  termination  will  be  approximately  the  same  generated  by  the  pump — i.e., 


Q = Q (3-7) 

o o 


Also  from  Eq.  (3-6),  if  t Fie  load  impedance  is  much  greater  than  the  source  impedance, 


Z 

o 


»ZS1  , then: 


(3-8) 


As  mentioned  previously,  a truly  anechoic  termination  can  only  be  approached  in 
liquid  systems.  Since  there  will  be  small  reflections,  the  easiest  way  to  measure  the 
incident  pressure  is  by  the  use  of  multiple  transducers.  The  recommended  stationary 
transducer  configuration  is  shown  in  Fig.  2-2.  With  the  use  of  the  OSU  Six  Transducer 
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Data  Reduction  Program  in  Appendix  C,  the  reflection  factor  of  the  termination  and 
the  incident  pressure  from  the  pump  can  be  determined.  The  maximum  allowable  re- 
flection from  the  termination  is  a reflection  factor  of  0.23  which  corresponds  to  a 
stadning  wave  ratio  of  4 dB. 

With  the  incident  pressure  from  the  pump  known  and  the  sonic  velocity  either  meas- 
ured or  estimated  by  the  computer  program,  the  pump  flow  ripple  can  be  calculated  using 
Eg.  (3-2). 

To  compare  the  magnitude  of  flow  ripple  for  one  pump  with  that  of  another,  the  flow 
ripple  of  the  first  and  second  harmonics  should  be  added  together  for  a particular  speed. 
If  this  is  done  for  several  speeds,  and  a linear  regression  performed  that  yields  a curve 
passing  through  the  origin,  then  a single  number  rating  can  be  used.  This  number  would 
be: 


Rating  = qjQy/N  (3-9) 

where:  qToT  = q1  + q2 

N = speed 

So,  in  order  to  calculate  the  approximate  flow  ripple  for  a certain  pump,  this  number 
would  only  have  to  be  multiplied  by  the  speed.  Obviously,  this  number  does  not  give 
the  purchaser  of  the  pump  enough  information  about  the  pump  to  adequately  predict  what 
it  will  do  in  a system;  however,  it  does  give  the  purchaser  a number  to  compare  with 
other  pumps.  A pump  with  a lower  number  would  be  one  that  would  generate  less  flow 
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ripple  and  vice-versa. 


The  test  procedure  for  determining  the  flow  ripple  of  a pump  is  given  in  Appendix  C. 


Pump  Impedance  Determination 

To  evaluate  the  impedance  of  a pump,  another  series  of  tests  must  be  performed  in 
which  a known  load  impedance  is  inserted  in  the  system  with  the  pump.  A simple  device 
that  can  be  used  for  a load  impedance  is  the  orifice. 

As  previous  tests  indicate  [9],  it  is  desirable  to  experimentally  determine  the 
characteristics  of  the  load  impedance  instead  of  just  assuming  it  a purely  resistive 
component.  To  evaluate  the  load  impedance,  the  test  procedure  outlined  in  Appendix  A 
should  be  followed.  In  this  case,  the  load  impedance  is  treated  as  a fluidborne  noise 
attenuator  with  the  transmission  loss  and  reflection  factor  being  determined  from  the 
data  taken  from  the  test. 

For  the  qeneral  case  of  a load  impedance,  Z-j.,  downstream  a distance,  L,  from  the 
pump  (See  Fig.  3-2.),  the  pressure  at  any  point,  d (from  the  load),  in  the  circuit  is 
given  by  [7]:  • 

P(d)  = Po6 (eY  d + p e d)  (3-10) 

. . -2Y L 

H-Pspte  j 

where*  P = pressure  that  would  occur  in  an  anechoicly  terminated 

o 

system 

p - termination  reflection  factor 
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t 


Since  P , P(d)  and  p are  either  measured  or  experimentally  determinable  quanti- 
o t 

ties,  p is  the  only  unknown  in  Eq.  (3-11)  or  Eq.  (3-12).  The  source  reflection 
coefficient,  p^,  is  related  to  the  source  impedance  by  the  followinq  equations: 

Z = Z (1  + P ) / (1  -p  ) (3-13) 

s o s s 

where:  Z = source  impedance 

Z = characteristic  line  impedance 
o 


The  source  impedance  has  both  a magnitude  and  phase  angle.  The  test  procedure 
for  determining  the  pump  impedance  is  given  in  Appendix  D.  Witn  both  the  source 
impedance  and  flow  ripple  determined,  the  fluidborne  noise  characteristics  of  a fluid 
power  pump  are  sufficiently  defined  to  predict  the  resultant  pressure  ripple  in  a system. 


TEST  SYSTEMS 

A schematic  of  the  anechoic  termination  test  system  was  shown  in  Fig.  3-1  . The 
anechoic  termination  was  constructed  out  of  41  meters  of  2.54  cm  0.  D.  steel  tubing 
and  39  meters  of  2.54  cm  0.  D.  steel  hose.  With  long  lengths  of  tubing  and  hose,  an 
anechoic  termination  is  accomplished  by  dissipating  the  pressure  ripple.  A valve,  to 
allow  pressurization  of  the  system,  and  a filter  are  located  downstream  of  the  anechoic 
termination.  Since  the  pressure  ripple  is  significantly  dissipated  by  the  time  it  reaches 
the  valve  and  filter  any  reflections  from  these  components  are  minimized.  Though  there 
is  dissipation  in  t fie  line,  it  is  assumed,  for  the  flow  rates  used,  that  the  effect  of  dissa- . 
pation  between  transducer  positions  is  negligible. 
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The  transducers  are  piezoelectric  pressure  transducers.  Each  transducer  is 
calibrated  to  within  ± 1%  accuracy.  The  transducers  are  mounted  in  individual  mounts 
so  that  the  tip  of  the  transducer  lies  flush  with  the  interior  wall  of  the  tubing.  The  trans- 
ducer spacing  is  the  same  as  shown  in  Fig.  2-3.  Because  of  the  generally  long  wave- 
lengths of  pressure  waves  in  hydraulic  oil,  the  transducers  must  be  spaced  far  enough 
apart  to  adequately  measure  the  standing  waves  in  the  system.  A more  thorough  dis- 
cussion of  the  use  of  stationary  pressure  transducers  for  fluidborne  noise  analysis  can 
be  found  in  Ref . [9  - 11]. 


Fig.  3-3.  "Close-Coupl  ing"  of 
Pump,  Pressure  T ransducer  & 
Orifice 


The  system  used  to  measure  the  transmission  loss  and  reflection  characteristics 
of  the  orifice  is  the  same  as  that  used  in  the  previous  chapter  (Fig.  2-1).  Measure- 
ments are  taken  both  upstream  and  downstream  of  the  orifice.  With  the  aid  of  the  OSU- 
Six  Transducer  Data  Reduction  Program  (See  Appendix  E.)  the  incident  and  transmitted 
pressure  can  be  determined,  as  well  as  the  reflection  factor  and  the  location  of  the 
standing  wave  maxima  and  minima  (v.hich  yields  the  reflection  factor  phase  angle)  up- 
stream of  the  orifice. 
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Shown  in  Fig.  3-3  is  the  orifice,  pump, and  pressure  transducer  arrangement  used 
to  evaluate  the  pump  impedance.  The  orifice  is  "close-coupled"  to  the  pump,  being 
less  than  1/20  of  a wavelength  away  from  the  pump. 

Test  Results 


The  procedures  outlined  in  this  report  for  measuring  the  fluidborne  noise  character- 
istics of  a hydraulic  pump  require  the  use  of  an  anechoic  termination.  This  section 
presents  the  results  of  evaluations  of  the  anechoic  termination,  the  flow  ripple  of  several 
pumps,  the  reflection  characteristics  of  an  orifice  and  the  internal  impedance  of  a pump 
for  several  frequencies. 


Anechoic  Termination 

Plotted  in  Fig.  3-4  is  the  reflection  factor  for  the  anechoic  termination  as  a function 
of  frequency.  Pressure  measurements  were  taken  with  six  stationary  pressure  trans- 
ducers. The  reflection  factor  was  determined  from  the  pressure  measurements  using 
the  OSU  Six  Transducer  Data  Reduction  Program. 
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Fig.  3-4.  Reflection  Factor  for 
Anechoic  Termination  Configuration 
Showing  the  Recommended  Maximum 
Acceptable  Reflection  Based  on 
4 dB  Standing  Wave 
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This  program  reconstructs  the  standing  wave  in  the  system  given  the  six  pressure 
measurements  and  frequency.  The  incident  pressure,  reflection  factor,  sonic  velocity, 
and  position  of  standing  wave  maxima  are  varied  until  a "best  fit"  combination  is  found 
for  the  data  input . 

The  largest  experimental  value  of  was  0.3878  at  100Hz.  This,  of  course, 
is  over  the  maximum  acceptable  value  for  the  reflection  factor.  This  was  true  of  several 
other  frequencies  as  shown  in  Fig.  3-4. 

As  the  flow  rate  and  frequency  increases,  the  general  trend  is  for  the  reflection 
factor  to  decrease.  The  frictional  dissipation  is  both  a function  of  the  flow  rate  and 
frequency.  As  either  the  flow  rate  or  frequency  increase  so  does  the  frictional  dissi- 
pation. Since  the  performance  of  this  type  of  anechoic  terrrination  depends  significantly 
upon  frictional  dissipation,  an  increase  in  this  factor  should  decrease  the  reflection  factor 
as  shown  in  Fiq.  3-4. 

Flow  Ripple  in  the  Line 

Three  identical  gear  pumps  were  used  to  obtain  an  average  value  of  experimental 
flow  ripple  in  the  line,  Q^,  from  a particular  model  pump  and  to  determine  a variation 
between  units  that  should  be  "identical."  The  flow  ripple  is  zero  at  zero  PPM,  and 
theoretically  should  increase  linearly  as  the  speed  of  a gear  pump  is  increased  [2], 

Based  upon  what  the  theoretical  values  should  be,  straight  lines  were  fitted  through 
the  origin  for  the  data.  Fig.  3-5  and  Fig.  3-6  show  the  data  and  the  fitted  lines  for 
the  first  two  harmonics  for  t Fie  three  pumps  used.  Because  of  the  high  reflection 
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Fig.  }-5.  Flow  Ripple  Versus  Pump 
Speed  Show i nr)  Experirnent.il  Values 
and  Best  Fit  Line  Through  Origin  for 
the  First  Harmonic 


Fig.  3-6.  Flow  Ripple  Versus  Pump 
Speed  Showing  Experimental  Values 
and  Best  Fit  Through  Origin  for  the 
Second  Harmonic 


encountered  at  the  low  flow  rates  generated  by  the  pumps,  several  data  points  for  both 
the  first  and  second  harmonics  were  not  used.  Since  the  termination  would  not  be 
considered  "anechoic"  at  those  particular  combinations  of  flow  rate  and  frequencies, 
the  experimental  values  of  flow  ripple  would  normally  be  erroneous. 

The  calculated  values  for  the  slopes  of  the  best  fit  lines  through  the  origin  are  listed 
in  Table  f-1  . A "t"  test  was  performed  on  the  data  to  test  the  hypothesis:  "the  fitted 

line  goes  through  the  origin."  The  results  of  the  test  are  also  shown  in  Tables  3—1  and 
t-2.  In  all  cases,  the  hypothesis  could  not  be  rejected  which  indicates  that  it  is 
acceptable  to  fit  this  set  of  data  with  a line  that  goes  through  the  origin. 
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The  fitted  line  for  a system  pressure  of  69  BAR  in  both  cases  had  a slope  larger 
than  the  ones  for  system  pressures  of  34.6  BAR.  However,  because  of  the  small  number 
of  data  points  in  each  case  and  the  larqe  variation  between  units,  it  was  not  possible  to 
distinguish  any  significant  difference  between  the  two  lines.  For  instance,  the  limit 
on  the  slope  of  the  second  harmonic,  34.5  EAR  line  was  2.65  x 10  ^ (8-  /min)/RPM  for 
the  lower  limit  and  1 .308  x 10~^  U/min)/RPM  for  the  upper  limit.  The  slope  of  the 
second  harmonic  69  BAR  line  was,  from  Table  3-1  , 8.52  x 10  k U /min)/RPM  which  is 
well  within  the  limits  of  the  aqcepted  slope. 


TABLE  3-1  . Slope  of  Flow  Ripple 
Line  and  Statistical  I nformation  for 
for  Two  Harmonics  and  Two  Pressures 


HARMONIC 

SYSTEM 

PRESSURE 

(BAR) 

SLqPE  OF 
LINE 

'V  Ml Nl  RPM 

CALCULATED 

t' 

1 

34  5 

1 48  « 10  3 

0 436 

690 

1 82  « 10  3 

0 602 

2 

34  5 

6 67  * 10  4 

0 707 

680 

8 52  « 10  4 

0 662 

FREQUENCY 

(Hz) 

REFLECTION 

FACTOR 

PHASE  ANGLE 
(DEGREES) 

200 

0 833 

1.49 

400 

0 9477 

5.78 

600 

0.9914 

17  39 

800 

0 9430 

6 68 

TABLE  3-2.  Reflection  Factor  and 
Phase  Angle  for  the  Orifice 
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These  prel  iminary  flow  ripple  tests  seem  to  indicate,  that,  for  a given  pump  model, 
there  can  be  a large  amount  in  variation  of  flow  ripple  between  units.  As  wide  a varia- 
tion as  4.  3 dB  (64.8%)  was  observed  in  the  900  RPM  second  harmonic,  69  BAR,  flow 
ripples  between  units  OSU-NP-42  and  OSU-NP-43.  The  largest  observed  standard 
deviation  for  any  one  data  point  was  less  than  0.6  dB  which  seems  to  indicate  that  the 
differences  between  the  various  pumping  units  was  due  to  something  inherent  in  the  pump- 
ing units  and  not  rreasurement  error. 

Pump  Impedance 

The  experimental  reflection  factor  and  phase  angle  of  the  orifice  are  listed  in 
Table  3-2.  The  phase  angle  of  the  orifice  was  calculated  by  determining  how  far  the 
standing  wave  maxima  was  shifted  away  from  the  orifice  toward  the  pump.  The  position 
of  the  standing  wave-  maxima  was  calculated  in  the  OSU  Six  Transducer  Data  Reduction 
Program  (See  Appendix  F.). 

The  experimental  reflection  factor  of  the  orifice  is  significantly  larger  than  that 
calculated  by  only  considering  the  resistive  component  of  the  impedance.  For  instance, 
at  a flow  rate  of  33.01  H/min.  and  pressure  drop  of  56.4  BAR,  the  expected  reflection 
factor  would  be  0.7421  if  only  resistance  were  considered. 

With  tine  reflection  factor  and  phase  angle  of  the  orifice  known,  the  orifice  was 
" lose -coupled"  to  the  pump  as  was  shown  in  Fig.  3-4.  The  pressure  was  measured 
, .irig  t‘  i transducer,  and  its  value  was  used  in  f qs.  (3-6)  and  (3-6)  to  determine  the 
i.urtt  r rfle<  lion  factor  and  impedance.  These  values  are  I isted  in  Table  3-3.  The 
hr  St  f.  jr  I ar  monies  were  usea  at  the  fastest  test  speed  (1200  RPM)  of  the  pump. 


1-12 


FUCOUCNCY 

(Hi) 

ft 

4>* 

z.** 

e* 

200 

0976 

14  110 

2 498 

84  4 

400 

0 601 

31  367 

0 829 

44  34 

600 

0 522 

37  177 

0 678 

40  90 

SOO 

0 702 

64  483 

0 447 

68  17 

TABLL  3-3.  Reflection  Factor 
and  Impedance  of  Pump  OSU-NP-43 
For  Various  Frequencies 


The  impedance  rf  the  pump  appears  to  decrease  with  increasinq  frequency  as  shown 
in  Fiq.  3-7  and  as  expected  with  an  analytical  model  of  a positive  displacement  pump  as 
proposed  by  Davidson  [Id).  In  the  analytical  model  the  pump  impedance  is  modeled 
proportional  to  the  inverse  of  the  frequency  squared  and  the  internal  volume  of  the  pump 
as  qiven  by  the  expression: 


where: 


Z 

s 


P c 


j w V 


P density  of  fluid  medium 

c sonic  velocity 

a)  frequency  in  rad/sec 

V internal  volume  of  pump 


(3-14) 


An  internal  volume  of  the  pump  OSU-NP-43  of  34.4  cm*  (2.1  in.'*)  and  the  averaqe 
sonic  velocity  of  1207. S m/sec  was  used  to  calculate  the  theoretical  impedance  of  the 
pump  sfiown  in  Fiq.  3-7.  The  larqest  observed  deviation  was  at  400H2  where  the 
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experimental  and  theoretical  values  of  impedance  were  off  by  less  than  3.5.  (This  is 
a dimensionless  number  because  the  logarithim  of  the  pump  impedance  divided  by  the 
characteristic  impedance  of  the  line  is  a dimensionless  quantity.) 


P 


. 

I 


. 
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Pump  Flow  Ripple 

With  the  impedance  of  the  pump  and  the  flow  ripple  in  the  anechoic  termination  known, 
it  was  then  possible  to  determine  the  pump  flow  ripple  using  Eq.  0-4).  This  was  done 
for  pump  CSU-NP-43  at  a system  pressure  of  69  BAR  for  two  different  anechoic  termina- 
tions. The  results  are  shown  in  Figs.  3-8  ana  3-9  for  the  first  and  second  harmonics, 
respective!  y . 

Again,  due  to  excessive  reflections  at  some  frequencies,  it  was  not  possible  to 
use  all  the  data.  That  is  why,  for  inst  ,nce,  at  the  first  harmonic,  data  is  not  plotted 
for  the  2.16  cm  I . D.  conduit  anechoic  termination  at  pump  speeds  of  600,  900,  and 
1 200  PPM . One  anechoic  termination  was  const  rue  ted  of  2.16  cm  I.  D.  steel  hose  and 
tubing,  while  the  other  was  constructed  of  1 .08  cm  I . D.  steel  hose  and  tubing. 

A linear  regression  was  performed  on  the  data  for  both  harmonics.  For  the  first 
harmonic,  the  correlation  coefficient,  r,  was  0.9789.  The  values  of  the  coefficients, 

A and  B,  for  the  regression  line,  A + B , were  A .410  H /min  and  B .0015  1/ min/ 
RPM.  For  the  second  harmonic,  r = 0.9485,  A = .400  t/min,  and  B .00031  fc/min/ 
RPM.  Because  of  the  large  correlation  coefficient,  it  would  appear  that  a linear  relation- 
ship exists  in  the  range  studied  between  the  free  volume  flow  ripple  (what  the  pump  actually 
generates)  and  the  speed  of  the  pump.  For  the  two  different  anechoic  terminations, 
that  data  at  any  one  point  for  the  free  volume  flow  ripple  was,  in  all  cases,  with. in  5% 
of  each  other.  This  would  indicate  the  validity  of  l q.  (3-4)  which  states  that  the  free 
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volume  flow  ripple  of  the  punv  is  independent  of  the  load  impedance.  Since,  in  the  data 
taken,  two  different  load  impedances,  were  used  will,  an  equivalent  free  volume  flow  ripple 
determined. 


hiq.  3-7.  Pump  I mpedance  as  a Function  of  Frequency 
(OSU-NP-43)  Nondimensional 
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I iq.  5-8.  Free  Volume  Flow  Ripple  as  a Function  of  Pump 
Speed — First  Harmonic 
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Fig.  3-9.  Free  Volume  Flow  Ripple  as  a Function  of  Pump 
Speed — First  Harmonic 
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CHAPTER  IV 

COMPUTER  PROGRAM 

The  purpose  of  this  chapter  is  to  discuss  the  theory  behind  the  OSU  Six  Trans- 
ducer Data  Reduction  Program  and  outline  how  this  program  can  be  used  to  assist  in 
analysis  of  the  fluidborne  noise  (FBN). 

Theory 

In  many  hydraulic  systems,  wave  propagation  in  lines  often  can  be  considered 
one  dimensional  ar.d  frictionless.  For  a simple  system  such  as  that  shown  in  Fig. 
there  is  a FBN  source,  such  as  a pump,  a transmission  line,  and  some  downstream 
impedance,  Zj,  that  has  a reflection  factor,  p 
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Fig.  4-1  . Schematic  of  Simple  Acoustic  System 
Consisting  of  Source,  Transmission  Line,  and  a 
Lead  Impedance,  Zj 


1-37 


The  pressure  ripple  that  the  FBN  source  generates  at  a given  frequency  can  be 


written  in  the  trigonometric  form: 


P = 

where:  P = 

o 

U)  = 

t 

k 

c = 

X = 


P Cos(wt  - kx) 
o 

amplitude  of  pressure  ripple 
angular  frequency  in  radians/sec 
time 

u 

c 

sonic  velocity 
distance 


(4-1) 


Lquation  (4-1)  describes  a pressure  wave  of  amplitude,  P , traveling  in  the  positive 
x-di  recti  on . 

When  the  pressure  wave  reaches  the  termination  impedance,  part  of  it  is 
back  toward  the  source,  and  part  of  it  will  continue  traveling  in  the  "positive 
The  portion  reflected  back  toward  the  source  can  be  described  as: 

p = p P Cos  (“t  + Kx)  (4-2) 

r t o 

Now  there  is  a "left"  and  "right"  traveling  wave  in  the  system  shown  in  Fig.  4-1  . 
The  wave  described  by  Eq.  (4-2)  will  travel  to  the  left  until  it  reaches  the  source,  at 
which  time  a portion  of  it  will  be  reflected  back  toward  the  right.  It  will  have  a magni- 
tude of  PT  P P . This  wave  will  then  travel  'own  the  pipe  until  it  is  reflected,  etc. 

T s o 


reflected 

x-direction." 
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The  sum  of  the  "right"  traveling  waves  incident  on  the  termination  is  given  by: 


P = P . Cos  ( cot  - kx) 
w 


(4-3) 


where: 


(4-4) 


P|  = Po(1  + PT  Ps  + ( pjp  + 1 PTPS)3  + * ' • ) 


Similarly  the  sum  of  the  "left"  traveling  wave  is: 


PL  = PR  COS  _ 


(4-5) 


where: 


(4-6) 

PR  = PTP|  = PTP0(1  + °T  ps  + (PTPS^2  + (pTps)j  + • • •> 

The  pressure  ripple  at  any  one  point  in  the  system  is  just  the  combination  of  the  "left" 
and  "right"  traveling  waves: 

P(x,t)  = PjCosfwt  - kx)  + PRCos  (w  t + kx)  (4-7) 

It  would  be  desirable  to  be  able  to  describe  the  magnitude  of  the  standing  wave 
in  the  system  as  just  a function  of  x.  In  other  words,  to  separate  Eq.  (4-7)  into  a 
magnitude  and  phase  component  of  the  form  [12]: 

P(x,t)  = P(x)  Cos  (dl  ) (4-8) 
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If  the  trigonometric  iaentities: 


and 


Cos  (wt  - kx)  = coswt  coskx  + sirwt  sinkx 


Cos  ( wt  + kx)  = coswt  coskx  - sinw  t sinkx 


are  substituted  into  Eq.  (4-7),  and  both  sides  squared,  t hie  resultant  reduces  to: 


P(x,t)  = [ P | ^ + P^2  + 2P|Pp  cos  2kx]'//2  cos  (i  I + ) (4-9) 

where: 

p = tan  ^ P I PR  tankx  (4-10) 

P + P 
I R 

Equation  (4-9)  is  ir  the  form  required  by  Eq.  (4-8).  The  magnitude  P(x)  of  the  standing 
wave  is: 


P(x)  = [P.?  + P 2 + 2P,P„  cos  2kx ] 1 /2  (4-11) 

I K IK 

P 

Since  pT  R_  this  substitution  can  be  made  into  Eq.  (4-11)  to  obtain: 

P, 

? 9 1/9 

P(x)  = P | [1  p -j.  + 2p-j-  cos  2kx]  ' (4-12) 


If  stationary  pressure  transducers  are  used  to  measure  theFBN  in  a system,  then  , 
at  least  three  measurement  locations  are  needed  to  solve  the  unknowns  ( P^  , p ^ and  x) 
in  Eq.  (4-12).  If  six  transducers  are  used  in  the  configuration  shown  in  Fig.  4-2, 


Tig.  4-2.  Six  Pressure  Transducers  Arbitrarily  Spaced 
in  a FBN  Measurement  System 

then,  the  pressure  ripple  measured  at  each  transducer  can  be  related  tc  x,  and  R by 
the  following  set  of  equations: 


P,  (x)  = P,2  [1  + P j2  + 2 PT  cos  2k  (x  -S})]W2  (4-13) 

p2<x)  = P(~  M + pj2  + 2 py  cos  2k  (x  -62)]'^^  (4-14) 

P3(x)  = P|“  [1  + py2  + 2 py  cos  2k  (x  (4-15) 

P4(x)  = P|2  [1  + Py2  + 2 Py  cos  2k  (x  -fi^)]1^2  (4-16) 

P5(x)  = P j 2 [1  + Py2  + 2 py  cos  2k  (x  (4-17) 

P&(x)  = P|2[1  + py"  + 2 py  cos  2k  (x  - |5f))l'/"’  (4-18) 

The  choice  of  lettering  transducer  4(b)  at  x was  completely  arbitrary.  The  value 


of  x gives  the  distance  that  transducer  4 is  relative  to  the  maxima  of  the  standing  wave. 
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Generally  speaking,  the  computer  program  finds  the  best  fit  of  Eg.  (4-15)  through 


Eg.  (4-18)  to  the  data  by  varying  x,  P-j. , Pj  , and  sonic  velocity. 

T he  Program 

The  following  are  needed  as  input  to  the  computer  program: 

(1)  the  measured  values  of  the  FBN  at  the  six  transducer  positions 

(2)  the  frequency  in  cycles/sec 

(3)  the  distances  6 ^ , 6 6 ^ , 6^,  and  6 ^ as  shown  in  Fig.  4-2. 

The  first  computational  step  in  the  program  is  converting  the  pressure  in  dB  to  PSI  . 
This  is  accon.pl  i shed  in  the  subroutine  PSI  . This  particular  subroutine  is  designed 
for  use  with  the  calibration  of  instrumentation  used  here  at  the  Fluid  Power  Research 
Center.  A certain  reading  is  taken  straight  from  the  charts  and  this  number  fed  as  input 
into  the  program.  This  number  needs  138.3  added  to  it  to  get  the  actual  pressure  in  dB . 
This  is  built  into  the  PSI  subroutine.  If  this  program  is  to  be  used  in  other  facilities, 
then,  this  subroutine  can  be  altered  to  suit  the  particular  needs  of  that  facility.  The 
pressure  in  PSI  is,  then,  calculated  by  the  subroutine. 

The  program,  then,  locates  the  measured  maximum  and  minimum  pressure  read- 
ings from  the  data  given.  These  are  calculated  for  two  reasons: 

(1)  It  gives  an  initial  estimated  value  of  the  reflection  factor,  p The  re- 

flection factor  is  related  to  the  actual  maximum  and  minimum  pressures  in 
the  pipe  at  a certain  frequency  by  equation: 


p 


(— Ei!£ 1)  / ( — max__  + ,} 

P p . 

min  rr.m 


(4-19) 


If  a certain  p and  p are  measured,  then,  the  actual  value  of  the  reflection 
max  mm 

factor  must,  at  least,  be  as  larqe  as  that  corresponding  to  the  measured  p 

max 

and  p . . 
min 

(2)  The  second  reason  p and  p . are  located  is  to  allow  an  initial  "ouessing" 

max  min 

point  for  the  value  of  x,  which  is  the  distance  that  transducer  4 is  from  the 
actual  maximum  of  the  standing  wave. 


As  has  been  stated  previously,  the  four  variables,  x,  , P ^ , and  the  sonic 
velocity  are  varied  in  order  to  best  fit  the  expected  curve  to  the  data.  The  best  fit 
is  found  by  "nesting"  four  do-loops  within  each  other  similar  to  that  shown  in  Fig.  4-3. 


SONIC 
VELOCITY 
00 -LOOP 


X- POSITION 
DO -LOOP 


INCIDENT  PRESSURE 
DO- LOOP 


REFLECTION  FACTOR 
DO -LOOP 


Fig.  4-  3 . Simple  F low  Diagram  of 
Important  Computational  "Nested" 
Do-Loops  in  Program 


f 
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A four-dirr.ensional  matrix,  in  essence,  is  searched  in  fixed  increments  until  a best  fit 


is  found.  The  incremental  size  step  on  the  sonic  velocity  is  45  m/sec  starting  at 
b45  m/sec  and  moving  up  to  1320  m/sec . The  position,  x,  is  varied  25  equal  incre- 
ments over  a wavelength. 


The  reflection  factor  can  be  varied  with  equal  increments  or  unequal  increments 
dependent  on  the  initial  value  of  the  estimated  reflection  factor.  For  a reflection 
factor  P yl.0.75  the  program  will  break  the  value  1 - P j into  twenty  equal  segments. 
If  p > 0.75,  the  program  will  use  a value: 


1 -P 


as  incremental  values  of  the  reflection  factor. 


The  incident  pressure  begins  with  an  estimated  value: 


P . = (P  max  * Pmin) 

4 


It  is,  then,  varied  as: 


(Pmax  + Pmin)  - P 


r + 1 


I ]*  0.4 
r 


The  use  of  the  computer  program  provides  an  approximation  of  the  parameters  needed 
to  define  the  standing  wave  in  a system.  It  is  only  as  good  as  the  theory  behind  it  and 
the  accuracy  of  the  pressure  transducers.  If  there  is  a significant  amount  of  dissipation 
in  a hydraulic  line,  then,  the  basic  equations  would  have  to  be  modified  to  include 
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dissipation.  in  that  case,  a new  unknown  would  be  entered  into  the  equations. 

However,  for  the  case  where  dissipation  is  not  a problem,  the  program  can  be  used 
with  only  minor  modifications  to  suit  a particular  facility  to  estimate  the  standing  wave 


parameter . 


CHAPTER  V 


CONCLUSIONS  AND  RECOMMENDATIONS 

Both  of  the  proposed  test  procedures  have  been  tested  in  the  laboratory  and  modi- 
fied to  incorporate  the  results  of  the  experimental  validation  as  well  as  the  inputs  from 
industry.  Eioth  procedures  form  valid  nuclei  for  industrially  acceptable  test  codes. 

The  test  code  for  evaluatinq  the  fluidborne  noise  generation  potential  of  hydraulic 
pumps  is  based  on  the  use  of  an  anechoic  termination  in  the  hydraulic  load  system. 

This  recommenaat ion  does  not  preclude  the  use  of  a procedure  which  requires  the  use 
of  a non-anechoic  termination  in  the  test  system. 

A test  code  under  consideration  by  the  NFPA  Project  Group  T3.9.24  incorporates 
two  non— anechoic  terminations  to  obtain  the  flow  ripple  and  impedance  of  a pump.  The 
code  was  developed  by  Davidson  [18].  Tests  are  currently  beiny  conducted  by  several 
individuals  on  the  NFPA  committee  to  test  the  repeatability  of  the  test  code.  The  pumps 
used  for  the  tests  were  provided  by  the  FPRC. 

It  is  recommended  that  the  test  code  developed  by  Davidson  for  the  NFPA  commi-tee 
be  compared  to  the  test  code  developed  by  FPRC  personnel.  The  basis  of  the  compari- 
son should  be  tests  per  formed  on  the  same  pumping  units.  Any  differences  in  either 
the  flow  ripple  or  impedance  for  the  two  codes  should,  then,  be  resolved. 


PA3s 


JVG7 
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Both  the  fluidborne  noise  attenuator  and  pump  fluidborne  noise  test  procedure  use 
the  digital  computer  in  data  reduction.  It  is  recommended  tf  it  the  computer  program 
continue  to  be  developed^*  pr  ovide  direct  < amputation  of  the  transmission  loss  of  FBN 
attenuators  and  FBN  characteristics  of  pumps. 

For  all  fluidborne  noise  measurements  and  predictions  in  hydraulic  systems,  a note 
of  caution  should  be  observed  regarding  the  magnitude  of  the  oscillatory  pressures, 
which  are  usually  significant  relative  to  the  mean  pressure  levels.  Nonlinearities  re- 
sulting because  of  the  high  noise/mean  pressure  ratios  will  become  more  evident  as 
we  gain  more  knowledge  in  this  area.  They  may  not  be  a major  concern. 

The  MERADCQM-OSU  Fluidborne  Noise  Attenuator  Test  Procedure  outl  ined  in 


this  report  should  be  submitteo  to  industry  for  review  . 
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APPENDIX  A 


Ml  THOD  FOR  EVALUATING  THE  PERI  ORMANCE 
OF  A FLUIDBORNE  NOISE  ATTENUATOR 


INTRODUCTION 

In  fluid  power  systems,  power  is  transmitted  and  controlled  through  a fluid  under 
pressure  within  an  enclosed  circuit.  Pumps  convert  mechanical  power  into  fluid  power. 
Sound  is  created  in  the  fluid  during  the  power  conversion  process.  In  fluid  power 
systems,  fluidborne  noise  can  be  controlled  with  the  use  of  fluidborne  noise  attenuators. 
The  attenuator  can  reflect  sound  energy  back  toward  the  source,  and  abscrb  sound  energy 
w it  hin  i t . 

The  need  fcr  a technique  to  evaluate  the  performance  of  fluidborne  noise  attenuators 
has  long  been  recognized  [12].  One  of  the  primary  considerations  in  the  selection  of 
a fluidborne  noif.e  attenuator  is  its  effectiveness  in  reducing  the  fluidborne  noise  in  a 
system . 

1 . Scope 

To  include  rating  basis,  test  methods  and  method  of  reporting  transmission 
loss,  pressure  drop,  and  impedance  of  fluidborne  noise  attenuators. 

2.  Purpose 

To  establish  a uniform  standard  for  measuring,  reporting  and  accurately  com- 
paring the  transmission  loss,  pressure  drop,  and  impedance  of  fluidborne  noise 
attenuators . 

3.  Terms  and  Definitions 

3.1  For  definition  of  Fluid  Power  terms  used  in  this  document , see  Ref.  [13]. 

3.2  For  definition  of  Acoustical  Terms  used  in  this  document,  see  Ref.  [14]. 

3.3  Anechoic  Termination.  A termination  in  a hydraulic  system  which  essen- 
tially reflects  none  of  the  incident  sound  energy  over  the  frequency  range 
of  interest . 

4.  Unit  of  Measurement 

The  International  System  of  Units  (SI)  is  used  in  accordance  with  Ref.  115] 
dB  relative  to  20  yPa. 

5.  Graphic  Symbols 

Graphic  symbols  used  are  in  accordance  with  Refs.  [16]  and  [17].  Where 
Refs.  [16]  and  [17]  are  not  in  agreement.  Ref.  [17]  governs. 


6. 


Test  Equipment 

6.1  Use  standard  measurement  devices  to  measure  oil  flow  , oil  pressure, 
pump  speed  and  oil  temperature. 

6.2  Use  one-third  octave  bund  or  narrow  band  analyzer  to  measure  fluid- 
borne  noise  levels. 

6.  1 Use  oscilloscope  to  measure  sonic  velocity. 

6.4  Use  oil  as  the  test  fluid. 


7. 


Test  Conditions  Accuracy 

7.1  Maintain  these  test  conditions  during  the  test  within  the  limits  specified 
in  T able  A-1  . 


TABLt  A-1 

Test  Conditions  Accuracy 


TEST  CONDITION 

CUSTOMARY 
U.  5.  UNIT 

S.  I.  UNIT 

MAINTAIN  WITHIN 
( * ) OF  ACTUAL 

FLOW 

USOPM 

H/MIN 

• 22 

PRESSURE 

PSI 

BAR 

- 22 

SPEED 

RPM 

REV  MIN 

t 22 

TEMPERATURE 

0 F 

° C 

18  t 2°C 

SONIC  VELOCITY 

FT/SEC 

M/SEC 

* 52 

8.  Test  Procedure 

8.1  Install  Anechoic  Termination 

8.2  Construct  Hydraulic  Circuit 

8.2.1  Install  in  the  circuit  oil  filters,  oil  coolers,  reservoirs,  and 
control  valve  as  needed  downstream  cf  the  anechoic  termination. 

8.2.2  Use  test  fluid  and  filtration  in  accordance  with  the  manufacturers' 
recommendation . 

6.2.  5 Install  inlet  and  discharge  line  diameters  in  accordance  with  the 
manufacturers'  recommended  practice.  Exercise  extra  care 
when  assembling  inlet  lines  to  prevent  air  leaking  into  circuit. 


I-S2 


8.  1 


8.4 


8.  8 


8.6 


Install  Pressure  Transducers 

8.3.1  I nstal  I two  sets  of  pressure  transducers  as  shown  in  Fig.  2-2 . 

8.3.2  Insure  that  transducer  spacing  is  the  sarre  as  in  Fig.  2-S  with 
spacing  tolerance  within  - 2%. 

8.3.  ] Place  a section  of  test  system,  equivalent  in  length  to  the  atten- 
uator , between  the  two  sets  of  pressure  transducers. 

Verify  Anechoic  Termination 

8.4.1  Set  system  pressure  to  maximum  test  pressure  for  FUN  attenuator. 

8.4.2  Set  speed  of  pump  so  that  the  fundamental  pumping  frequency  is 
at  100  Hz. 

8.4.3  Record  pressure  transducer  measurenr ents  as  a function  of  fre- 
quency in  dB,  for  first  three  harmonics  of  fundamental  pumping 
frequency . 

8.4.4  Measure  a "tare"  pressure  drop  vs.  flow  rate. 

8.4.8  Measure  sonic  velocity  per  OSU-FPRC  procedure  [12]. 

8.4.6  Repeat  8.4.2,  8.4.3,  8.4.4,  and  8.4.8  for  fundamental  pumping 
frequencies  of  128,  160,  26V,  and  315Hz. 

8.4.7  Insure  that  the  range  of  dB  readings  ^.4  dB  for  all  test  frequencies. 

8.4.8  Record  ali  data  in  data  table  similar  to  Table  A-2. 

Measure  Attenuator  Performance 

8.8.1  Remove  the  section  of  test  system  equivalent  in  length  to  the  atten- 
uator between  the  sets  of  pressure  transducers. 

8.8.2  Install  FBN  attenuator. 

8.8.3  Repeat  step  8.4.1  . 

8.8. 'i  Repeat  step  8.4.2. 

8.8.8  Repeat  measurements  in  8.4.3. 

8.8.6  Measure  flow  and  pressure  drop  across  attenuator . 

8.8.7  Relieve  system  pressure  to  a minirrum  pressure. 

8.8.8  Increase  system  pressure  to  that  in  8.8.3. 

8.8.9  Repeat  measurements  in  8.8.8  and  8.8.6. 

8.8.10  Repeat  steps  8.8.7  and  8.8.8. 

8.8.11  Repeat  measurements  in  8.8.9. 

8.8.12  Repeat  steps  8. 5 . 4 up  to  and  inc luding  8. 8 . 1 1 for  fundamental 
pumping  frequencies  of  128,  160,  267,  and  318  Hz  as  indicated 
in  T able  A-2  . 

Invert  attenuator  and  repeat  steps  8.8.3  through  8.8.12. 


Calculations 

9.1  Obtain  pressure  drop  versus  flow  due  to  the  attenuator  by  correcting 
8.5.6  with  8.4.4. 

9.2  Calculate  the  transmission  loss  and  impedance  of  the  attenuator  from  data 
taken  in  8.8. 


Data  Presentation 

10.1  Report  pressure  drop  versus  flow  rate  in  a form  similar  to  that  in  Fig.  A-1. 

10.2  Report  attenuator  transmission  loss  as  a function  of  frequency  as  shown  in 
Fig.  A-2. 

10.3  Report  attenuator  input  reflection  factor  as  a function  of  frequency  as  shown 
in  F ig.  A-3 . 

10.4  Report  attenuator  output  reflection  factor  as  a function  of  frequency  as 
shown  in  Fig.  A-3. 


TABLE  A-2.  Data  Table  for  FBN  Attenuator  Test 


Date: 


Attenuator  I 


run  p Speed: 
Fluid: 


System  Pressure: 
Fluid  Viscosity: 


Pump  I D. 
Flow  Rate: 


Tare  Pressure  Drop: 


Pressure  Diff. 


Fundamental 
Pumping  Frequency: 


I »c 


20  JO  40  W 60 


FLOWf  GPM) 


Fig.  A-1.  Pressure  Drop; 
Versus  Flow  Rate 


TRANSMISSION  .OSS  148) 

VS 

fRlOuf NCY  (HZ) 

OSO-fA-iT 
Mi*  -L-Z104 
»MMP-OSU-NP  4 

wr»c 

SONIC  V€lOOTV.  i?§S  M/S 

THCO*«T<Al 

• f *P»RIM»N’Al 


'RfOCWNCv  'Ml 


Fig.  A-2.  Transmission  Loss 
Versus  Frequency 


t 


O 


Ntri  ACTION  FACTOR 

»• 

FREQUENCY 
OSU  i T 
THEORETICAL 

I ROM  OSU  SI*  TRANSOUCf  R REDUCTION  PROGRAM 
DEDUCED  FROM  TL  Eflin  91 


MIL  L ?>04 
PUMP  OSU  NP  4 
T • 27* 

SO NC  VELOCITY  <?«  M/S 


too  no  900  960  400  4 SO 

FREQUENCY  (MZ) 


Fig.  A-3.  Reflection  Factor 
as  a T unction  of  Frequency 


11.  Summary  of  Designated  Information 

The  following  designated  information  is  needed  when  applying  this  recommended 
standard  to  a particular  application  or  use. 

11.1  Transmission  loss  spectrum 

11.2  Pressure  drop  versus  flow 

11.3  Input  reflection  factor  spectrum 

11.4  Output  reflection  factor  spectrum 
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APPENDIX  B 


PUMP  FLUIDBORNE  NO  I SE  CHARAC TER  I S T I CS  TEST  PROCEDURE 

As  discussed  in  Chapter  III,  knowledge  of  two  pump  characteristics  is  necessary 
in  order  to  adequately  define  the  pump's  acoustical  performance  in  any  given  system. 
They  are: 

1 . The  flow  ripple  output  of  the  pump  (analogous  to  the  current  output  of  an 
electrical  current  source). 

2.  The  pump's  internal  impedance  (analogous  to  the  internal  impedance  of  an 
electrical  current  source). 

The  test  procedure  for  evaluating  the  pump  flow  ripple  output  is  given  in  Appendix 
C.  Appendix  C gives  the  test  procedure  for  evaluating  the  pump's  internal  impedance. 
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APPENDIX  C 


PUMP  FLOW  RIPPLE  EVALUATION  TEST  PROCEDURE 


1 . Mount  pump  and  pressure  transducers  with  anechoic  termination. 

2.  Test  anechoic  termination: 

a.  Set  mean  system  pressure  to  maximum  rating. 

b.  Insure  that  the  range  of  dB  readings  £.4  dB  at  each  frequency  over 
the  desired  frequency  range. 

3.  Measurement  of  pump  flow  ripple: 

a.  After  establishing  condition  2(b),  test  the  pump  at  the  following  speeds 

and  maximum  pressure:  one-fourth,  one-half,  three-fourths , and 

maximum  speed.  Record  the  first  three  harmonics  for  each  setting. 

b.  Repeat  3(a)  for  one-half  maximum  pressure. 

4.  Reducing  and  reporting  the  data: 

a.  From  the  tests  in  Part  2,  use  the  OSU  Six  Transducer  Data  Reduction 
Program  to  determine  the  reflection  factor  for  the  anechoic  termination. 
The  data  should  be  plotted  in  the  same  format  as  Fig.  A-3. 

b.  From  the  tests  in  Part  3,  use  the  OSU  Six  Transducer  Data  Reduction 
Program  to  determine  the  pressure  ripple  in  the  pipe  and  the  sonic  veloc- 
ity in  the  fluid. 

c.  Determine  the  value  of  the  pump  flow  ripple,  , using  Eq.  (3-2)  and 
Eq.  (3-4).*  The  flow  ripple  should  be  plotted  as  a function  of  the  speed 
of  the  pump  with  two  different  harmonics.  The  overall  value  of  the  flow 
ripple  can  be  obtained  for  each  speed  by  adding  q^  and  • A value  of 

* See  Appendix  D. 


PFECSDITC  PAGE  ELANK-NOT  f ILMaD 
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q/N  should  be  obtained  by  performing  a linear  regression  on  the  plot  of 
q-j-Q-j-  vs.  speed. 


APPENDIX  D 


PUMP  IMPEDANCE  EVALUATION  TEST  PROCEDURE 

1 . Mount  pump  and  pressure  transducers  with  ar.echoic  termination. 

2.  Test  anechoic  termination: 

a.  Set  mean  system  pressure  to  maximum  rating. 

b.  Insure  that  the  range  of  dB  readings  <.4  dB  at  each  frequency  over 
the  desired  frequency  range. 

3.  Measurement  of  incident  pressure,  P : 

o 

a.  After  establishing  condition  2(b),  test  the  pump  with  the  anechoic 

termination  at  maximum  pressure  and  one  of  the  speeds  used  in  the  Flow 
Ripple  Evaluation  per  Appendix  C.  Record  the  first  three  harmonics 
from  each  transducer. 

4.  Obtain  an  orifice  which  will  produce  a pressure  drop  equal  to  the  maximum 

pressure  rating  at  the  pump  flow  corresponding  to  the  speed  chosen  in  3(a). 

Install  orifice  in  test  system  per  Fig.  2-3. 

5.  Test  orifice  as  an  attenuator  per  Appendix  A.  Use  pump  speed  per  3(a). 

Reduce  data  per  Appendix  A and  obtain  p(  for  the  first  three  harmonics. 

6.  "Close-couple"  pump,  transducer,  and  orifice  per  Fig.  3-4. 

7.  Test  system  at  speed  per  3(a).  Record  first  three  harmonics  from  trans- 
ducer (P  (d)). 

8.  Reducing  and  reporting  t fie  data: 

a.  From  the  results  of  Part  2,  use  the  OSU  Six  Transducer  Data  Reduction 
Program  to  determine  the  reflection  factor  (p  (()  for  the  anechoic  termin- 
ation at  the  first  three  harmonics.  Record  on  Table  D-1  . 
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b.  From  the  results  of  Part  3,  use  the  OSU  Six  Transducer  Data  Reduction 
Program  to  determine  the  incident  sound  pressure  from  the  pump,  P , 
for  the  first  three  harmonics.  Record  on  Table  D-l  . 

c.  Reduce  data  from  Part  5 per  Appendix  A to  determine  p for  the  ori- 

to 

fice  for  the  first  three  harmonics.  Record  on  Table  D-l  . 

d.  Record  P(d)  data  from  Part  7 on  Table  D-l  . 

e.  Using  Eq.  (3-6)  determine  the  pump  reflection  factor,  p . Record 
on  T able  D-l  . 

f.  Using  Eq.  (3-7)  determine  the  pump  impedance,  Z . Record  on  Table 

D-1  . 


I 

f 


APPE  ME  I ^ ' 

1 HI  OSU  SIX  TRANSDUCE  P DA  I A Rt  DUCT  It 'N  PROGRAM 


The  OSU  Six  Transducer  Data  Reduction  Program  is  1 computer  pi'  gram  w»  a ' 
uses  the  sound  pressure  data  supplied  by  the  station. u y prc  s ir  > tin  ■ Ju<  s to  p» 
v ide  a best  estimate  of  th<  ire  distribul  t<  I th«  /sten 

accoustical  characteristics. 

Inputs  to  the  program  are  the  sound  pressure  level  readings  (in  dP  re  1 1 y 

from  the  analyzer  charts  and  the  pumping  frequency  (fundam.  nta!  nr  harmonn  i under 

investigation.  From  this  information  the  proqran  dedui  es  the  mi  .dent  . • md  pres-  ■ i 

(I  ) , t 1 ' reflec  t ion  factor  ( P ) , the  estimated  »oni<  velo<  ity  (<  t , the  po  it  ion  of  tl  t 
i s 

standing  wave  maximum  (X),  the  error  squared  between  measured  pre azures  and  "best 
fit"  pressures,  and  the  "best  fit"  press  ures  at  the  transducer  locations. 

The  program  is  listed  on  pp.b6-  76.  Page  77  is  an  example  print-out  of  the  pro- 
gram's output  for  two  frequencies. 
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O o o 


1 RELEASE.  2.0 


Ml  NP 


DAT  E 


/63  J5 


12/55/53 


C 

C 

C 

L 


SUBROUTINE  M INP( P , PM  I N,UM IN  ,0 ) 


♦ ♦ 
» THIS  SUBROUTINE  LOCATES  The  MEASURED  minimum  of  THE  ♦ 

♦ STANDING  mAVE  in  the  conduit  ♦ 

* P M I N = M I N I MU  M MEASURED  PRESSUPE  * 


DIMENSION  PI  10)  ,01 10) 

IF (PI  1 ) .LE  .PI  2 ) ) GO  TO  1 
PMI N=P|2 I 
DMIN=D( 2 ) 

GO  TO  2 
1 PMIN  = P l 1 ) 

DM  I N=D I 1 ) 

2 IFIPM IN.LE .PI  3)  ) GO  TO  3 
PMI N=P<  3 ) 

DM  I N=D I 3) 

i 1EIPMIN.LE.PI4I)  GO  TIJ  -♦ 
PM! N=P(4 ) 

UMIN= J.O 

4 IFIPMIN.LE.PI5))  GO  TO  5 
P m l N=  P 1 5 ) 

DM 1 N=-Ul 4) 

5 IF IPM IN.L  E .P (6  ) ) GO  TO  6 
PMI N = P I 6 I 

DM  I N=-D I 5) 
t>  RETURN 
E NJ 


C 

c 

c 

c 

c 

c 


SUBROUTINE  XMA XI C , S , YM , X 1 ,F R , S ON ) 

HmmiKIHtttHHIHHHHUtHtHHIHHIItHHHHttttHHt 


THIS  SUBROUTINE  CALCULATES  THE  VALUE  OF  X,  THE  POSITION 
OF  THE  MAXIMUM  PRESSURc 


OIMtNSION  Cl  10)  ,S(  101 
REAL  Kl 

Wll=YM*IC(5l*Cl4l-2) 

■ 12=  2*CI 31  -C I 2)-  Ctl) 

■ 21=  YM*(S(5)*S(4)  ) 

■22=  S(2)*S(1)-2*S<3) 
■=( Wl 1 ♦ ■ 12 I/IM21 tW22 ) 
Kl= (2*3. 1415 9*FR /SON) 
X1=I ATANI M) ) /( 2*K1 I 
RETURN 

END 
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Gl  RELEASf  2.0 


PR  EU 


DATE  - 763  3 6 


12/55/63 


C 

C 

c 

c 

c 

c 

c 

c 

c 


2 

3 


1 


A 


20 

21 


SUbROUT I NE  PREU( R E . PR . S N , X N , PNR , FK, Dl 

DIMENSION  XN  ( 10) ,RE(101,PR(10l,SN(10l,PNR(lJ),D(10I.X(10» 


♦ THIS  SUBROUTINE  CALCULATES  THE  BEST  FIT  PRESSURES  AT  THE  ♦ 

♦ SIX  TRANSDUCER  LOCATION:,  TO  BE  PRINTED  JN  THE  OUTPUT  OF  ♦ 

♦ THE  PRUGRAM  ♦ 

♦ t 
♦♦444444444444444444444fr44+444f44444444444444444444444f444444444*4 


l 3-0  .0 

IFISNI1I.EC.231  GO  TO  A 
R = RE  I 1 1 *RE  (1)4-1 
Ol=A»3.1Al59  *FR/SN( l) 
XIII = COS (01* IXNI1 1 -C ( l III 
X( 21 =COS<0 l* l XNl  11-0(21) 1 
x ( 3)=C0$ 1U1* (XNl 1)—  0(3))) 
X ( 5) =CGS (01*  (XN( U *D(A ) ) 1 
X(  A)=COS  <D 1“ ( XN( 1 ) ) ) 

X (6 ) = CUS (01  * (XN( 1) *D( 5 ) I I 

OG  i 1=1,6 

PN= 2 * X ( 1 ) «RE ( 1 J 4« 

IFIPN.LT.Z3)  GO  T 1 ? 

uj  TO  3 

PN=- PN 

CUNT  I NUE 

PNl  =Suk  T ( PM 

PNR ( I ) = P R ( 1 ) «PN1 

C °NT  1 NUE 

60  TO  21 

CC'NT  | NUE 

00  20  J=  1 ,6 

PNR ( J) =0." 

CONI INUE 
CUNT INUE 
KE  TURN 
tND 
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I RELEASE  2.  0 


PCAl 


U A IE 


/63  3 S 


12/SS/B3 


C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 


3 

10 


SUBROUTINE  PCAL(REF.X,FNS.PC) 
DIMENSION  XI  10), PCI  10) 


♦ 

♦ THIS  SUBROUTINE  CALCULATES  THE  ESTIMATED  VALUE  OF  THE  ♦ 

♦ PRESSURE  AT  THE  Six  TRANSDUCER  LOCATIONS  GIVEN  THE  ESTI-  ♦ 

♦ MATEO  VALUES  OF  INCIDENT  PRESSURE,  POSITION  OF  THE  ♦ 

♦ STANDING  WAVE  MAXIMA,  SONIC  VELOCITY,  AND  THE  REFLECTION  ♦ 

♦ FACTOR.  ♦ 

♦ ♦ 


Y l = REF*REF  +1 
DO  10  1=1,6 
Y2=2*REF  *Xl  I ) 

Y 3 = Y 1 ♦ Y 2 

IF! Y3.LE.0)  GO  TO  2 

YA=SURT I Y 3 I 

GO  TO  3 

Y3=- Y3 

YA=SORT(  Y5  ) 

CONT INUE 
PCI  I )=PNS*YA 
CuNI INUE 
RETURN 
I- NO 


j 

i 


I 


c 

c 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 


20 


SUBROUTINE  CMPRE «PC,/2.VRl) 
JI ME  NSI ON  PC (10)  , 22 ( 10  ) 


♦♦♦♦♦♦*< 

♦ 


♦ 

this  subroutine  compares  the  values  of  the  pressure  ♦ 

FROM  THE  BEST  FIT  OPTIMIZATION  WITH  THE  MEASURED  ♦ 

VALUES  AND  DETERMINES  THE  SUM  OF  THE  VARIANCES  FOR  THAT  ♦ 
PARTICULAR  FIT.  THIS  WILL  THEN  l,E  LUMPARtU  TO  THE  2 S ♦ 

PEP  CENT  ERROR  VARIANCE  IN  SUBROUTIN  REFP  OR  IN  THE 
MAIN  PROGRAM,  CEPENCING  UPON  wHcPE  IT  WAS  CALLED  FROM  ♦ 


Y2-0.0 
DO  20  K=1 ,6 
Y*22(R  )-PC(K  I 
Y1 * V **2 
V2=  VI ♦ V2 
CUNT INUE 
VR1 * Y2 
RE  TURN 
ENO 
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GI  RELEASE.  2.0  REGV  OATE  = To 335 

XNI  41=XNI  3 ) 

XNI 3 )=XN( 2) 

X N I 2 > =X  1 

G 3 TO  6 

2  IF IVRI.GE.UI 3) > GO  TO  3 
Ul  5)  = U<4) 

Gl 4) =01 3) 

U I 3 ) = VRl 
PRIS)  = PR (A  I 
PR|4 J=PRI  1) 

PR  I 3 I = PN  S 
RE  I !>  1 =RE  I 4 1 
REI4  ) = PE<  3) 

R E I 3 ) =RF  2 
SMI  5 ) = SN I 4 ) 

SM(4)=SNI 3) 

SNI 3 ) = SON 
XNI  3 ) = XNI 4 I 
X Ml  4 ) = XN I 3 ) 
xMl  3 I = X 1 
GO  TO  6 

3 IKVR1.GE.UI4I  ) GO  TO  4 
Ul  5)  =UI  4) 

U14»=VRI 
PR|*>)  = PR  ( 4 I 
PR  I *♦)  =PNS 
Rf I 5 I =k£ 14  I 
Rt 14) = RF  2 
SNI 51= SNI 4) 

S N 1 4 ) = SON 
XNI 5) = XNI 4) 

AN  I 4 1=  XI 
GO  TC  6 

4 IK  VRl.GE.Ul  5)  I OC  TO  5 
Jl 5 ) = VR  1 

PR  I 5 I = PMS 
RE(5)=RF2 
SNI 5 I = SON 
XNI  5 I =X1 
!>  CONTINUE 
6 CONTINUE 
RETURN 
END 


12/33/5  3 
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r>  0 r. 


R E L c A S E ? . 0 


RE  GV 


DATE 


It )>3*> 


12/SS/S  3 


C 

C 

C 

C 

C 

C 

C 


S'Jb-i  JUT  | NE  REGVI VR l ,U .PNS. PR.RF2 ,RE . SN.SON.xl  , XN  I 

♦ ♦ 

♦ This  IS  THE  SUBROUTINE  »Hfc  R E THE  BEST  FIT  VALUES  CF  l NC  I - ♦ 

♦ JENT  PRESSURE  IP  IN ) . REFLECTION  FACT  MSONI.  ERROR  SU  JAR  ED  ♦ 

♦ CVR1)  . AND  POSITICNIXil  ARE  STORED.  uNLV  THE  BEST  FIVE  ♦ 

♦ VALUES  ARE  STORED.  ♦ 


D I ME  NS  I ON  U(10I,PRU0I,KEI1D),SNU0I 

DIMENSION  XNIIOI 

IF IVR1 .GE.UI 1 I I GO  TO  1 

U1 SI =U14I 

U I 4 I =U(  3) 

U l 3 ) =U I 2 ) 

Ul c) =UC 1 I 

um  = vRi 

PR  IS  I =PR  <4  I 
PRI A) = PR  13) 

PR ( 3 l = PR (2  I 
PR (2  I =PR I 1 ) 

PRI  1 MPNS 
R E I S I = R E l 4 J 
Rt!4)=REl3> 

RE(  3) =R  E l 2 1 
F E I 2 I *RE I 1 1 
RE  U I = RF 2 
SNI S)=SNl4) 

SNlA )=SN( 31 
SMil  * SN 1 2 I 
SN I 2 1 = SN  I 1) 

SNI1  i = SON 
XNC  SI = XN ( A I 
XN  I A ) = XN  I 3 ) 

XN(3)=XN(2) 

XNC  2 ) >XN( 1 I 
XN  C l I = X I 
GO  T 0 6 

I IFCVR1. GE.UI  2)  I GO  TG  2 
U l S )=U( A ) 

U I A I -Ul 3) 

Ul  3 1 =Ul  2 ) 

U 12 l-VR  I 
PRISI  =PR I A I 
PRI A)=PRI 31 
PRI 3 J*PR (2  I 
PR  I 2 I * PN  S 
REIS)=RE(A) 

RE  I A ) *RE 1 3 ) 

RE(3I*REI2I 
REI 2 ) = RF  2 
SN(S)-SNIA» 

SNI  A I * SNI  3 ) 

SNI 3 )*SNI 2 ) 

SNI 2 I *SON 
XNI  5)  * XN  I 41 
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RELEASE  2.0 


COR 


DATE 


loi  3S 


12/SS/S3 


C 

C 

L 

C 

C 

C 

C 

i. 

c 


SuoHOUTINt  CORIX.FR, SCN.D.Xl) 


♦ ♦ 

♦ THIS  SUBROUTINE  CALCULATES  THE  VALUtS  OF  THE  CCSINES  OF  ♦ 

♦ TRANSDUCER  LOCATIONS  P-  EL  AT  1 V t TO  THE  POSITION  OF  THE  ♦ 

♦ STANDING  wAVe  MAXMA  XI, ONCE  XI  HAS  SEEN  ESTIMATED.  F 

♦ v 


DIMENSION  XI  10)  ,U<  10) 
01=4*3. l4l59*FR/S0N 
x( II = COSI D1 * ( Xl-OI 1 I I I 
Al 21=  COSI D l * I Xl-DI 211  ) 

X I 3 I = C0S(U1*1X1-DI 311  I 

XI  4)  = CO  SC  01  *X  1 ) 

XI  5>  )=  COSI  D 1*1 XI *01 41)) 
X16 ) = CUS 1D1  ♦ (X 1 »D(S ) ) I 
RE  TURN 
ENO 


,1  RfltASE.  Z.O 


KEEP 


DATE  = 7o  33S  1 2/SS/S3 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 


30 


fj 


4-t 


3 


31 

7 

SO 

SS 


10 


SUBROUTINE  kEFP<P|N,RE,PNS,RER,VAKU,Xl,.),FR,SON.  Z2, PM Ax. PMIN, PR, 
iU iSRi xn  I 

0 1 ME  NS  ION  X 110  1 ,0( 101  , RE  ( 10  ) , PS  l 1 0 1 , l 2 l 1 0) , PC  1 10 J 
DIMENSION  XN110) 

U 1 ME  NS  I ON  Pk  l 10)  .U<  10)  , >N(  IC» 


RITHIN  THIS  SUBROUTINE,  THE  INCIDENT  PRESSURE  AND 
THE  REFLECTION  FAcOR  ARE  VARIED  , ONE  AT  A TIME,  IN 
OROER  TO  ARRIVE  AT  THt  VALUES  OF  THE  St  KESPtCTEU 
PARAMETtRS  THAT  JEST  PIT  THE  DATA. 


tHttttl 


CALL  CUR  ( X.FR, SON.D, X 1 I 
N MAX =10 

DO  10  J=1,NMAX 
CALL  PCAL 1REF, X, PNS.PC 1 
CALL  ..MPRE  (PC,  12  ,VR1  I 
RF  2 = R L F 

IF (VR1.L  E.VARU)  GU  TO  30 

GO  T 0 t) 

CALL  KEG  VI VR 1, U, PNS.PR  ,«F2 , RE , SN, SON, XI ,X, ) 
NMX= 1 0 
ib= 0. 7" 

R I N= ( 1-RFF ) /NM  X 

DO  SO  K= 1 , NMX 

I F ( R F 2 • L T.ZS)  GU  TO  4^ 

■<F2=(  l*KF2)/2 

oC  TC  3 

«FZ  = RtF<-K»RI N 

CALL  PCAL (RF 2. X.PNS.PCI 

CALL  CMPRE (PC, Z2 ,VR2 I 

IF ( VR2.LE . VR  1)  GO  TO  4 

GO  TO  SS 

VKl = VR2 

IF( VR1.LE.VAR0I  GO  TC  31 
GO  TO  7 

CALL  RE GV ( VR 1 ,U , PNS , PR  , RF2 , RE , S N, SON , X 1 , xN  I 
CUNT  I NUE 
CCNT INUE 
CONT I NUE 

V l« ( ( PMA  X4PMIN1 Z2-PNS1 *0.40 

PNS  = PNS ♦ V 1 

CONTINUE 

RETURN 

t ND 
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O r> 


t 

C 

C 

C 

C 

C 

c 


SUBROUTINE  MAXP(P,PMAX,  JHAX.J) 


. + ♦ 

♦ THIS  SUBROUTINE  LOCATES  THE  HEASUREO  MAXIMUM  OF  THE  ♦ 

♦ S T A MO  I NO  IN  THE  CCNUUIT  ♦ 

♦ PMAX=MAX 1 MUM  MEASUREO  PRESSURE  ♦ 

♦ ♦ 

DIMENSION  PI  101 .01  10) 

IFIP111.GE.P12))  GC  TO  1 

PMAX=P12) 

l)MAX  = 01  2 ) 

GO  TO  2 
1 PMAX=P<  1 ) 

UMAX  = 01  1 ) 

2 IE  I PMAX . GE .P  13 ) I GC  TO  i 
P MA  X - P ( 3) 

0MAX=D( 3 ) 

3 1 FIPMAX.GE.PIA  ) I GO  TO 
PMAX=  P I A ) 

uMAX=0.0 

A I E I PMAX.GE . P ( 5 I I GU  T C 5 
PMA  X=  P(  5 ) 

UMAX  =-U( A I 

5 IF (PMAx.GE.P(b) ) GO  TC  o 
PMAX=  P ( 6 ) 

UMAX =-015) 

G RETURN 
rNE) 


C 

c 

c 

c 

c 

c 

c 

c 

G 


SUB ROUT I Nt  SUMOl 22 .VARU) 
01  ME. MS  ION  22(10) 


♦ ♦ 

♦ This  SUBROUTINE  CALCULATES  THE  SUM  of  THE  VARIANCES  IE  ♦ 

» THERE  MERE  A 25  PER  ufcNT  ERROR  IN  THE  FIT  PE  A CUMvt  TO  ♦ 

♦ THE  DATA.  THIS  GIVES  THE  MAXIMUM  ALLOWABLE  ERROR  IN  * 

♦ ATTEMPTING  TO  FIT  A CURVE  TU  THE  DATA  ♦ 

♦ VAKO  = SUM  OF  THE  VARIANCES  EUR  A 2s  PER  CENT  cRRI.R  ♦ 


w2  = 0 .0 

UO  1 1=1  ,b 

W = 22< I ) *0.  1 
Ml = M**2 
M2  =M  1 ♦ M2 
CUNT INUE 
V AR  U = m2  *6 .25 
RE  TURN 
END 
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1 RELEASE.  2.C 


FUNC 


DATE 


Tb  i ib 


12/bb/bi 


C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 


SUB ROUT  I me  F UNCI J 1 ,FR , SQN.U ,C , S J 
DIMENSION  C 1 101  , CUO)  , S(  10  ) 


♦ this  subroutine  calculates  THE  VALUES  of  the  sines  anu  ♦ 

♦ COSINES  FOR  THE  DISTANCE  A PARTICULAR  TRANSDUCER  IS  ♦ 

♦ FKCM  TRANSDUCER  A.  ♦ 

♦ 0< it  .. . .D( b> = DISTANCES  T1  , T2 , Ti , Tb . AMO  T6  ARE  FROM  T A ♦ 

♦ F R=  F RE  OU  EN  CY  ♦ 

♦ SQN  = SON I C VELOCITY  ♦ 

♦ ♦ 


Ol  = 2*<?*FK/SUN)*DlJll»3.1<,lbS 
ul J1  ) = CQS(D1  ) 

S(  J1  ) =SI  N(  01  ) 

RETURN 
c NO 


C 

c 

c 

c 

c 

L 

c 

c 

c 

c 

c 

c 

c 


SUbROUTI NF  P S I < P . J . P S , 12 ) 
DIMENSION  PI  IP), PS ( 10)  , L2 1 101 


+ 

THIS  SUBROUTINE  CCNVERTS  THE  PRESSURE  REAUINU  FROM  THE  ♦ 
TRANSDUCERS,  READ  IN  DEuluELS  CN  I HE  CHART,  TO  PS  I . ♦ 

THE  FOLLOWING  PARAMETERS  ARE  USE  ) ♦ 

P(JI»THE  PRESSURE  IN  DECIBELS  ♦ 

PS( J Is  THE  PRESSURE  IN  PSI  ♦ 


miHiH 


2 = <Pl JI-A2. FI/20 

l 1= 10»*Z 

22  I J J -2 .9*21 

PSI J)=22< JI*Z2I J) 

RETURN 

E ND 
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.1  RELEASE  .2  .0 


MA  IN 


,)MI 


loi)!) 


12/bb/bi 


I 


C 

c 

c. 

c 

c 

c 

c 

L 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

(. 

c 

l 

c 

c 

L 

c 


1 


6 


2 


60 


A 


u I M»  NS  I UN  u ( L 0 ) ,HE(10),Z2(IJ)>PSI  1 0 ) , P I 1 0 i « C ( 10)  tSI  10) 
U I ME  N S l ON  U(  1 ■>)  , SNt  1 0)  , P»<  1 0) 

o i me ns  ion  xi  n i , xni  io) , pnri  ioi 


♦ ♦ 

♦ THIS  PKCCKAM  ESTIMATES  THE  REFLECTION  E AC  T TP  t INCIDENT  * 

*■  PRESSURE,  POSITION  OF  THE  STANDING  WAVE  MAXIMA,  AND  THE  ♦ 

♦ SONIC  VELOCITY  FOR  THE  STANDING  WAV!.  IT  IS  DESIGNED  TO  bE  * 

♦ USED  WITH  DATA  FPCM  j I X STATIONARY  PRESSURE  TP  AN  SDUCER  S . ♦ 

♦ THE  NUMBER  OE  SETS  OE  DATA  MOST  »E  READ  IN,  PLUS  ThE  DIS-  + 

♦ TANCtS  BETWEEN  TRANSUUCtRS  1 AND  A,  2 AND  A,  3 AND  A,  S * 

♦ AND  A,  AND  6 AND  A.  THt  PRESSURE  DATA  FRuM  I HE  T PANS-  ♦ 

♦ DUCERS  IS  THEN  READ  IN  AS  wElL  AS  THE  TREUUENCY  FOP  cALh  * 

F LAT  ASET.  ♦ 


THE  MAJOR  PARAMETERS  U St U ARE: 

Dll)- DISTANCE  FROM  TRANSDUCER  l T ] A 
012) 'DISTANCE  Fk CM  TRANSDUCER  2 TO  A 
01  3) 'U I STANCE  From  TRANSDUCER  3 T ^ A 
UIAl'OISTANCE  from  transducer  A TJ  5 
01 6) 'Dl STANCE  FPCM  TRANSDUCER  A TO  b 
Xl'RELAT IVE  POSITION  UF  TRANSUUU  p a 
STANDING  WAVE  MAXIMA 
RE  F =RE  FLfcC  T I ON  FACTCR 
SON' SON  1C  VFLOCITYIIN  M/StL) 
PMAX'MEAjUkF J STANDING  WAVE  MAXIMA 
PM  I N 'Mb  A SURE  D STANDING  WAVE  MINIMA 
PNS= INCIDENT  PRESSURE 

U-FRROR  SwUARECl VAR  IANCE I FUR  ESTIMAI 


TRIM  I HE 


ED  CURVE 


F IT 


RFAUlS.l)  N 
FURMAT 112) 
wR  I T E I 6 , 3 3 ) 

33  FORMAT  I 1 H 1 ) 

H f A J I B , 6 I I 0 I I 1 ) , I 1 ' 1 , 5 ) 
FuRMA  T< SIF8. 2 ,2X) ) 

00  3 N 1 ' 1 , N 

RE  ADI S ,2 ) CPI  I I, l'l*b),FR 
FORMA  I I/IF  8. 2.2X I ) 

DJ  60  M 3 = 1 . 6 

XNIM3)  =0.0 

HE  I M3) =0. 0 

PK<  M3)=0  .0 

U!  M3) =1000000 

SNI  M 3) '0. 0 

CUNT INUE 

00  A J=1 ,6 

CALL  PS1 (P.J.PS.Z2) 

CONT INUE 

CALL  MAXPI 22  ,PMAX, CMAX , 0) 
CALL  M|NPIZ2,PMIN,UM1N,0) 
SWR'PMAX/PMIN 
CALL  SUMDl Z2 ,V ARO) 

DU  10  Ml=l, 16 
S0N'6a6*MI*A5 
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RtLEASE  2.0  MAIN  DATE  - 76333  12/33/53 

DO  70  N5=  1 . 2 5 

*EF=IS*R-1 1/  ISMR*1  ) 

PNS  = ( PMA  X*  PM  | N I / A 
fiAL  = S0N/(50«FR) 

XI =0M| N» 1 N5-1 I *HAL 

CALL  REEP(PIN,RE,PNS.REE,VARD.XI,E),ER,SCN,ZZ,  PM  AX,  PM  l N,  PR,U, SN, AN ) 
70  CONTINUE 
10  CONTINUE 

wR  l TE  I u.  3M) 

35  FORMAT  I5X , // //  ) 

*P(TE<6.377  E« 

37  FORMAT l l HO  , • FREUUf  NC  Y I S • , 2 X ,F  o , J , 2X  , • HZ  • , /) 

»RITE(6»36)  VAkJ 

3fc  FORMAT1  1H0,'  «AXI  MU'1  A L L rR  A o IE  ERROR  IS  *.iX,F8.2J 
v.ALL  PKEJIPF.PR,  SN.XN.PNR.f  K.DI 
*RI TC 16,22)  I Z2 I Ibl, 16=1,6) 

22  FORMAT  I 5 X , ’MEASURE  U PP E S SUR E S • , 3 X ,6 1 F d . 2 , 3 A ) I 

<*R  I T E |6 , 21  ) IPNR  l I 5) , I 5= 1, 6) 

21  FORMAT (3X, ‘REST  FIT  PR  ES  SUR  ES  ■ . 3 X , 6 I Fb  . 2 , 3X  ) ) 

«R 1 TE I 6. dO) 

dO  FORMAT (IX, 'REFLECTION  FACTOR • , 3X , • I NC I Ut NT  P R t S3 oP E • , 2 X , • E RF  Oh  SOU 
»AREU'  ,2X  SONIC  VELOCI  TV  ,dX,  ' POSIT  ION  I’E  X‘,/l 
JO  63  M2  = l , 5 

RRI  T E l o ,20  ) RE  IM2  ) ,PR  (M2  ) ,U  (M2  I , SN  ( M2  ) , XN  I M2  ) 

20  F IRMA! ( 5>X,Fo.4 , l2X,Fo. 2,dX,Fd.2  ,6X ,Ed.2  ,3X  ,Fo  .2,  / I 
65  CONTINUE 
3 CONTINUE 
NlUP 
END 
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